1/2 


AD-A161  899  CLUSTERING  AND  ORDERING  IN  III-V  ALLOVS<U>  WASHINGTON 
UN IV  ST  LOUIS  MO  SEMICONDUCTOR  RESEARCH  LAB 
C  M  WOLFE  ET  AL  21  JUL  85  UU/SRL-59582A-9 
UNCLASSIFIED  AFOSR-TR-85-0931  AFOSR-82-0231  F/G  11/6  NL 


1.0 


l.l 


|m 

III2'2 

■  36 

Ih 

1 2.0 

I  1.8 

4_ 

1.6 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  1963-A 


,  \  ,N  , 

.V'V*' 


.**  w’Vi  I 


WASHINGTON 
UNIVERSITY 
IN  ST  LOUIS 


*F<XR-TR  8  5  -  0  f 


CLUSTERING  AND  ORDERING  IN  Hl-V 


SEMICONDUCTOR  RESEARCH  LABORATORY 
Washington  University 
Saint  Louis,  Missouri  63130 


31  July  1985 


Dt;c 

kSLLC  ?>;*• 

1 3  wsf 

'• f  , 


m  a 

“4 


1  June  1984  to  31  May  1985 

ANNUAL  SCIENTIFIC  REPORT  NO.  WU/SRL-59583A-9 

Air  Force  Office  of  Scientific  Research 

Building  410 
Bolling  Air  Force  Base,  DC20332 


one  File  copy 

Grant  No.  AF0SR-82-0231 


*  •  "  "'-lbllo  reler.cn  ; 
-t  itievu  oaiiuited. 


The  United  States  Government  is  authorized  to  reproduce  and  distribute  this  report 


for  Governmental  purposes. 


85  11  12  138 


Unclassified 


7; 

V 

VC 


9 


i\ 


f ; 


fJ. 


c-' 

t 


i 


[:• 
t  •' 

1 


l 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  _  ,  _  ^  2.  GOVT  ACCESSION  NO. 

1EOSRTR.  8  5-  0  931  , 

3.  RECIPIENT'S  CATALOG  NUMBER 

09? 

«.  TITLE  (end  Submit) 

Clustering  and  Ordering  in  IIl-v  Alloys 

5.  TYPE  OF  REPORT  h  PERIOD  COVEREO 

Annual 

1  June  84-31  May  85 

6.  PERFORMING  ORG.  REPORT  NUMBER 

WU/SRL  59583A-9 

7.  AU  T HOR( s) 

C.M.  Wolfe,  M.W.  Muller,  S.  Julie  Hsieh, 

3.  CONTRACT  OR  GRANT  NUMBERfi; 

AFOSR-82-0231 

9.  PERFORMING  ORGANIZATION  name  and  aodress 

Washington  University 

Box  1127 

St.  Louis,  MO  63130 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

2306/BI 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Air  Force  Office  of  Scientific  Researh 
Building  410 

Bollina  AFB ,  DC  20332 

12.  REPORT  DATE 

31  July  85 

13.  NUMBER  of  pages 

133 

u.  MONITORING  AGENCY  name  ft  AODRESS  (if  different  from  Controlling  Of/ice; 

IS.  SECURITY  CLASS,  (a!  chit  report) 

Unclassified 

15a.  DECLASSIFICATION.  DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (oi  this  Report) 

The  United  States  Government  is  authorized  to  reproduce  and 
distribute  this  report  for  Governmental  purposes.  A 

.  „«Mlo  releaflM 

Approved  for  p '  lted. 

distribution  unll»i 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  ( Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

ZnSnP2,  InxGai_xP,  AfxGal-xAs'  GaAs,  hetero junctions ,  interfacial 
energy  gaps,  generalized  Wannier  functions,  alloy  formation, 
long-range  ordering,  primitive  layered  structure. 

20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identity  by  block  number ) 

-“Staggered  lineup  hetero  junctions  are  expected  to  exhibit  tun¬ 
neling  assisted  optical  transitions  across  the  interfacial  energy 
gap  Ej.  In  the  staggered  lineup ,>Ej  is  smaller  than  either  of  the 
energy  gaps  of  the  constituent  semiconductors.  In  this  work,  we 
examined  two  staggered  lineup  hetero junctions,  ZnSnP^/GaAs  and 
Ini_xGaxP  GaAs .  Below  bandgap  emission  and  absorption  were  ob¬ 
served  in  i.nese  structures,  allowing  the  experimental  determination  „ 

DD  1  j*N  7J  1473  COITION  OF  t  MOV  «5  IS  OBSOLETE 


ITnrrl  agg-i  f  i  oH 


V 


. 


I  SECURITY  CLASSIFICATION  OF  THIS  PACE  fWhwi  Date  Entered) 


of  l£f.  Good  agreement  was  seen  between  this  value  of  and 
that  from  the  predicted  band  lineup. 


Advances  in  technology  have  made  possible  the  fabrication 
of  rapidly  varying  heterostructures  which  hold  the  promise  of 
important  applications.  We  develop  a  set  of  approximate  treatments 
of  electron  states  in  a  variety  of  layered  heterostructures.  The 
approximations  are  all  based  on  the  concept  of  one-band  generalized 
Wannier  functions .  Following  a  discuggimmf  Llie-"ra'Mdity  pf  ihia)  ■  ■■ 
representation,  we  apply  it  to  an  evaluation  of  the  bound  states 
in  a  narrow  quantum  well  in  GaAs,  which  clearly  demonstrates  the 
mixing  of  main  and  satellite  valley  states  as  well  as  the  contribu¬ 
tion  of  evanescent  states.  We  also  evaluate  the  states  of  a  super¬ 
lattice  in  a  model  structure  of  up  to  20  quantum  wells.  As  a  final 
example  we  discuss  the  application  of  generalized  Wannier  functions 
to  the  matching  of  electronic  states  at  a  hetero junction  between  two 
model  band  structures  with  different  effective  masses,  and  compare 
the  formalism  with  alternative  approaches  to  this  problem. 

Recent  III-V  alloy  formation  models  indicate  that  negative 
charge  transfer  energy  can  overcome  positive  bond  distortion  energy 
to  stabilize  long-range  order  or  compound  formation.  Although  a 
Export  of  a  layered  ordering  in  AlxGai_xAs  tends  to  confirm  this 
result,  we  have  as  yet  obtained  no  convincing  evidence  for  such 
ordering  in  InxGai_xP,  Compound  formation  in  these  alloys  may  be 
critically  dependent 'em  growth  temperature  and  kinetic  factors. 


Section 


Page 


TABLE  OF  CONTENTS 


1.  Introduction .  1 

2.  Staggered  Lineup  Hetero junctions .  2 

2.1  ZnSnP2/GaAs .  13 

2.2  In1_xGaxP/GaAs .  60 

2.3  C-V  Intercept  Method .  70 

2.4  Conclusions .  77 

2.5  References .  81 

3.  Spatially  Varying  Energy  Bands .  83 

3.1  One-Dimensional  Wannier  Model .  86 

3.2  General  Wannier  Model .  93 

3.3  Simple  Examples .  97 

3.4  "True"  Heterostructures . 106 

3.5  Conclusions . 112 

3 . 6  References . 115 

4.  Long  Range  Ordering . 119 

4.1  Recent  Alloy  Formation  Models . 119 

4.2  X-ray  Measurements . 121 

4.3  Layered  Structure . 125 

4.4  Conclusions . 129 

4.5  References . 129 

5.  Personnel . 130 

6.  Publications . 131 

7.  Talks . 13  2 

.'.-’7  Tfnnr-  QyjTQ  -  -■  ,  -„TC  i 

TROTS*  1  :  • 

Thi  i  tec-  -  •  v 

appro-.  ••  *o»*  tnhi  t«  •  •  '  ; 

IMstr-  :»:>i  •  .,U 

'  ■■  TTh  ,;t 


.  Technical  Inforaat ion  Division 


INTRODUCTION 


Although  there  is  substantial  interest  in  III-V 
semiconductor  alloys  for  electronic  and  optoelectronic 
devices,  relatively  little  is  known  about  the  distribution 
of  atoms  in  these  materials  or  the  effects  of  non-random 
distributions  on  device  performance.  In  many  of  these  alloys, 
however,  non-random  distributions  of  atoms  are  expected  to  be 
important  from  thermodynamic  considerations.  This  is  par¬ 
ticularly  true  at  the  low  temperatures  commonly  used  for 
epitaxial  growth.  The  available  evidence  suggests  that 
short-range  clustering  of  like  atoms  or  short-range  ordering 
of  unlike  atoms  can  produce  device  problems  such  as  excess 
noise  and  leakage  current,  premature  voltage  breakdown,  and 
lower  carrier  mobilities.  Long-range  ordering  of  unlike  atoms, 
however,  could  potentially  yield  III-V  ternary  compounds  with 
properties  superior  to  their  parent  alloys.  In  epitaxial 
growth,  substrate  effects  such  as  lattice  match  and  non¬ 
equivalent  sublattice  sites  are  expected  to  have  a  strong 
influence  on  these  phenomena.  Also,  two  of  the  important 
III-V  alloys,  GaxIn]__xp  a^d  GaxInl-xAs'  can  be  epitaxially 
grown  matched  to  GaAs  and  InP,  respectively,  at  compositions 
near  the  point  of  greatest  ordering  probability.  The  objective 
of  this  work  is  to  investigate  various  aspects  of  clustering 
and  ordering  in  III-V  alloys,  including  the  interfacial 
properties  of  the  materials. 


-„v„\ 
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2.  STAGGERED  LINEUP  HETERO JUNCTIONS 

There  are  three  flat  band  hetero junction  lineups 

with  non-zero  band  offsets.  These  flat  band  lineups  are 
shown  in  Figure  1,  where  the  three  types  of  lineup,  the 
straddling,  the  staggered,  and  the  broken-gap  lineups 
are  named  according  to  terminology  from  a  previous 
publication  [1].  The  lineups  shown  are  called  flat  band 
lineups  to  indicate  there  is  no  band  bending  resulting 
from  deviations  from  local  charge  neutrality.  Band  bending 
in  the  equilibrium  and  forward-biased  lineups  wi". 1  be 
considered  later.  The  straddling  lineup,  which 
is  the  most  common,  results  when  the  valence  band  maximum 
of  the  larger  bandgap  material  lies  below  that  of  the 
smaller  bandgap  material.  The  less  common  staggered  lineup 
has  both  band  edges  of  one  material  below  the  respective 
band  edges  of  the  other  material.  The  fairly  rare  broken- 
gap  lineup  is  that  having  the  conduction  band  minimum 
of  one  semiconductor  below  the  valence  band  maxima  of 
the  other. 

The  band  lineup  we  will  be  concerned  with  in  this 
work  is  the  staggered  lineup.  The  reason  this  lineup 
is  so  interesting  is  because  of  the  energy  gap  that 
exists  between  the  lower  conduction  band  edge  and  the 
higher  valence  band  edge  at  the  hetero-interface.  This 


interfacial  energy  gap,  Ej,  as  seen  in  the  middle  of 


Figure  1,  is  smaller  than  the  bandgaps  of  either  of  the 
constituent  semiconductors.  This  is  the  essential  feature 
of  the  staggered  lineup  he tero junction  that  has  led  to 
discussions  of  this  lineup  being  a  potential  source  of 
below  bandgap  emission  and  absorption  [2,3]. 

Several  different  lattice-matched  materials  systems 

have  been  studied  to  determine  the  viability  of  the 

staggered  lineup  hetero junction  as  a  below  bandgap  emitter 

or  detector.  E.J.  Caine  and  co-workers  have  investigated 

p-Al0  4aIn0  52As/n-InP  grown  by  molecular  beam  epitaxy 

(MBE)  [4].  They  observed  below  bandgap  electroluminescence 

at  1.4K  which  peaked  at  0.97eV  for  a  pulsed  current  density 
2 

of  4.5  A/cm  .  They  used  this  peak  energy  value  as  roughly 
equal  to  Ex,  and  found  reasonable  agreement  with  their 
predicted  value  of  Ex=l.leV. 

Another  materials  system  which  has  been  studied  to 
investigate  the  optical  properties  of  staggered  lineup 
hetero junctions  is  p-In-  Ga  P/n-GaAs  grown  by  liquid 

X  "  X  X 

phase  epitaxy  ( LPE ) .  S.  Julie  Hsieh  and  co-workers  reported 
seeing  a  below  bandgap  peak  in  the  photovoltaic  response 
at  room  temperature  [3].  The  half -power  point  of  the 
peak  was  at  1.31eV  and  corresponded  closely  to  their 
predictions  for  the  interfacial  energy  gap.  The  predicted 
band  lineup  for  this  hetero junction  was  actually  a 
straddling  lineup  rather  than  a  staggered  one. 


Experimentally  these  lineups  were  measured  to  be  staggered. 
This  may  be  the  result  of  the  particular  orientation  of  the 
GaAs  substrate  used,  {111},  which  may  have  interfacial 
dipoles  that  shift  the  lineup  from  a  straddling  to  a 
staggered  one.  For  reasons  explained  in  a  later  section, 
we  also  explored  these  p-In1-xGaxP/n-GaAs  hetero junctions 
as  potentially  superior  below  gap  emitters. 

The  third  materials  system  that  has  been  examined 
to  study  the  staggered  lineup  hetero junction  is 
p-2nSnP2/n-GaAs.  At  room  temperature  we  observed  below 
bandgap  peaks  in  both  the  emission  and  absorption  spectra 
which  we  believe  may  be  due  to  the  below  gap  mechanisms 
discussed  in  the  next  section.  Experimental 
values  of  EI  obtained  from  these  data  agreed 
well  with  the  theoretical  prediction  of  Ej  obtained  from 
the  calculated  band  lineup.  The  ZnSnP2/GaAs  hetero junction 
was  the  one  investigated  in  the  greatest  depth  and  so 
these  results  will  constitute  the  major  portion  of  this 
work. 

To  understand  some  of  the  possible  mechanisms  for 
below  bandgap  absorption  and  emission  in  the  staggered 
lineup  hetero junction,  let  us  consider  the  appropriate 
band  diagrams  for  p-ZnSnP2/n-GaAs .  The  flat  band  lineup 
for  this  hetero junction  was  calculated  using  results  from 


capacitance-voltage  (C-V),  Hall,  photoconductivity,  and 
photovoltaic  measurements ,  which  will  be  discussed  in 
detail  in  Section  2.1.  Figure  2  reveals  that  the  flat 
band  lineup  for  ZnSnP2/GaAs  is  staggered  with  the 
conduction  band  minimum  of  the  p-type  ZnSnP2  below  that 
of  the  n-GaAs.  This  configuration  has  been  referred  to 
as  the  p  (down)/n  (up)  band  lineup  [2].  Since  ZnSnP2  is 
grown  unintentionally  p-type,  this  configuration  is  the 
only  one  we  have  been  able  to  achieve  with  this  material. 

To  investigate  the  opposite  configuration,  the 
n  (down)/p  (up)  lineup,  we  examined  the  p-In.  Ga  P/ n-GaAs 

A 

hetero junction.  These  results  are  discussed  in  Section  2.2 

Now  examine  the  equilibrium  band  lineup  of 
p-ZnSnP 2 /n-GaAs  in  Figure  3.  Charge  has  transferred  to 
equilibrate  the  electrochemical  potential  throughout  the 
heterojunction.  This  charge  transfer  creates  depletion 
regions  just  as  in  the  p-n  homo junction  case,  except  now 
the  band  offsets  create  barriers  and  form  potential  wells 
on  adjacent  sides  of  the  interface.  This  diagram  is  drawn 
in  exaggerated  proportions  to  aid  in  the  analysis  of  this 
section. 

Consider  the  effect  of  illumination  on  this  hetero¬ 
junction  with  zero  applied  voltage.  A  potential  mechanism 
for  below  gap  absorption  might  involve  the  excitation 
of  an  electron  out  of  a  state  in  the  valence  band  well 


of  the  n-GaAs  to  a  virtual  state  in  the  gap.  Since  the 
electron  considered  quantum  mechanically  is  a  wavepacket, 
part  of  the  wavepacket  of  an  electron  in  a  virtual  state 
near  the  interface  may  overlap  unfilled  states  in  the 
conduction  band  well  of  ZnSnP2.  In  other  words,  the 
electron  may  tunnel  from  the  virtual  state  into  unoccupied 
states  in  the  ZnSnP2  well.  The  energy  of  the  photon 
absorbed  therefore  corresponds  to  the  energy  difference 
between  the  initial  valence  band  state  and  the  virtual 
state  in  the  gap  of  GaAs .  This  energy,  as  can  be  seen 
from  the  diagram,  is  less  than  the  energy  gaps  of  either 
of  the  constituent  semiconductors.  The  below  gap  absorption 
is  then  detected  as  the  electrons  in  the  conduction  band 
well  and  holes  created  in  the  valence  band  well  tunnel 
through  their  respective  barriers  and  drift  to  contribute 
to  a  negative  generation  current. 

A  second  mechanism  would  involve  the  tunneling  of 
an  electron  out  of  the  valence  band  well  to  a  virtual 
state  followed  by  excitation  to  the  conduction  band  through 
photon  absorption.  This  second  mechanism  is  equivalent 
to  hole  tunneling  and  so  is  much  less  probable  than  electron 
tunneling  due  to  the  larger  hole  effective  mass.  This 
is  because  the  tunneling  probability  of  a  carrier  through 
a  barrier  decreases  roughly  exponentially  with  effective 
mass.  For  example,  the  tunneling  probability  for  an 


electron  tunneling  through  a  triangular  barrier  is, 

*  |  va 

P(fi)  =  exp  (~4(2mJ  A-£ - ]  ,  (1.1) 

3qhE 

where  E  is  the  constant  electric  field,  £  is  the  electron 

energy,  as  is  the  energy  barrier  height  above  £  ,  and 
* 

in  is  the  electron  effective  mass. 

To  examine  a  potential  mechanism  for  below  bandgap 
emission  look  at  the  hetero junction  in  forward  bias  in 
Figure  4.  The  applied  voltage  is  indicated  in  terms  of 
the  potential  energy  difference,  qV,  between  the  quasi-Fermi 
levels  in  the  field  free  regions  of  the  two  materials. 

In  forward  bias  minority  carriers  diffuse  over  the  dimin¬ 
ished  barriers  and  are  temporarily  confined  in  the  potential 
wells.  The  electrons  confined  in  the  conduction  band 
well  of  the  p-ZnSnP2  could  then  tunnel  out  to  virtual 
states  in  the  gap  of  GaAs  and  then  radiatively  recombine 
with  a  hole  confined  in  the  valence  band  well  of  the 
n-GaAs.  Another  mechanism  would  involve  the  electron 
emitting  a  photon  and  dropping  from  the  conduction  band 
well  to  a  virtual  state  in  the  gap  of  ZnSnP2.  The  electron 
could  then  tunnel  to  recombine  with  a  hole  confined  in  the 
GaAs  valence  band  well.  This  second  mechanism  involves 
hole  tunneling  and,  as  is  the  case  for  absorption,  would 
be  a  less  likely  mechanism.  A  previous  publication  [2] 
suggested  that  the  opposite  staggering  arrangement,  the 


n(down) /p(up)  lineup,  might  prove  to  be  the  most  strongly 
emitting  lineup  since  in  this  case  accumulated  majority 
carriers  instead  of  injected  minority  carriers  are  involved 
in  the  radiative  recombination.  Since  we  have  not  been 
able  to  grow  n-type  ZnSnP2,  we  obtained  this  lineup  using 
p-In. _  Ga  P/n-GaAs  he ter o junctions.  The  results  of 
emission  experiments  on  these  devices  will  be  examined 
in  Section  2.2.3. 
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2.1  ZnSnP2/GaAs 

2.1.1  Material  Properties  and  Device  Construction 

Before  considering  the  construction  of  the  devices 
used  in  these  investigations,  it  will  be  useful  to  discuss 
some  important  material  properties  of  the  somewhat 
unfamiliar  semiconductor  ZnSnP2.  ZnSnP2  is  a  II-IV-V2 
ternary  compound  which  can  have  either  the  disordered 
sphalerite  or  the  ordered  chalcopyrite  crystal  structure. 
The  sphalerite  structure  for  ZnSnP2  is  essentially  the 
cubic  zinc-blende  structure  having  a  face  centered 
cubic  (fee)  sublattice  of  P  atoms  and  a  random  fee 
sublattice  of  the  Zn  and  Sn  cations.  The  chalcopyrite 
structure  is  tetragonal  because  the  cation  sites  are 
nonequivalent  and  ordered.  The  chalcopyrite  structure 
is  shown  in  Figure  5,  where  the  tetragonal  axis  is  desig¬ 
nated  as  the  c  axis.  In  ZnSn?2,  the  chalcopyrite  structure 
shows  no  tetragonal  distortion,  unlike  some  chalcopyrites , 
so  that  c  =  2a.  Using  double-crystal  x-ray  diffraction 
measurements  the  lattice  constant  of  ZnSnP2  was  determined 
as  5.6507  A,  which  results  in  a  lattice  mismatch  with 
GaAs  of  less  than  0.44%.  The  energy  bandgap  of  ZnSnP2 
was  determined  to  be  1.45eV  through  differential  analysis 
of  photoconductivity  data  on  ZnSnP2  layers  grown  on 
semi-insulating  (SI)  GaAs  substrates. 
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The  device  structure  used  in  this  work,  shown  in 

Figure  6,  consisted  of  p+-ZnSnP2  layers  grown  on  (211}.n+- 

GaAs  substrates.  The  ZnSnP2  layers  were  grown  by  LPE  using 

an  open  tube  system.  These  layers  were  grown  by  Gary  A. 

Davis  and  the  details  of  the  growth  were  reported 

earlier  [5].  The  ZnSnP2  was  grown  unintentionally  p-type. 

19  -3 

The  hole  concentration  was  determined  to  be  p  =  5  x  10  cm 

from  Hall  measurements  performed  on  layers  grown  on  SI 

GaAs  substrates.  The  electron  concentration  of  the  GaAs 

substrate  was  determined  using  Hall  measurements  to  be 
18  -3 

n  =  10  cm  .  The  bandgap  of  the  GaAs  substrates  was 
observed  by  photovoltaic  measurements  to  be  1.42eV.  We 
also  observed  evidence  of  Zn  diffusion  from  the  epitaxial 
ZnSnP2  into  the  GaAs  substrates.  Delineation  etches  of 
ZnSnP2  layers  grown  on  SI  GaAs  revealed  both  a  growth 
interface  and  a  diffusion  front  at  about  1  um  into  the 
substrate.  This  indicated  that  the  Zn  diffusion  was 
converting  a  thin  region  of  the  SI  GaAs  substrates  near 
the  interface  to  p-type.  When  the  ZnSnP2  was  grown  on 
n+-GaAs,  the  Zn  diffusion  was  not  sufficient  to  convert 
the  GaAs  to  p-type.  In  this  case  only  a  single  growth 
line  was  visible  at  the  interface  after  the  delineation 
etch.  Microprobe  analysis  indicated  outdif fusion  of  Ga 
into  less  than  1 ym  of  the  ZnSnP2  but  with  no  type 
conversion.  The  ZnSnP2  layers  used  in  these  devices  were 
fairly  thin,  about  3um.  The  diodes  were  defined  using 
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standard  photolithographic  techniques  and  were  etched 

-3  2 

down  to  an  area  of  2  x  10  cm  using  Br-methanol.  Finally, 
10  mil  Sn  spheres  were  alloyed  to  both  materials  since 
Sn  forms  ohmic  contacts  with  both  p-ZnSnP2  and  n-GaAs. 

The  GaAs  substrates  used  in  this  work  had  the  {211} 
orientation.  This  choice  of  substrate  orientation  was 
particularly  important  for  two  reasons.  Firstly,  previous 
work  [6]  revealed  that  the  {211}  orientation  of  the  GaAs 
suppresses  all  but  one  orientation  of  the  c-axis  of  the 
chalcopyrite  ZnSnP2  layer.  The  c-axis  of  the  ZnSnP2  can 
have  three  different  orientations  with  respect  to  the 
cubic  axes  of  the  sphalerite  GaAs  substrate.  Where  regions 
having  different  c-axis  orientations  meet,  an  antiphase 
domain  boundary  is  produced.  Antiphase  domain  boundaries 
contribute  to  strain  and  scattering  and  so  degrade  the 
mechanical  and  electrical  properties  of  the  material. 

Back  reflection  Laue  photographs  of  ZnSnP2  layers  grown 
on  several  different  substrate  orientations  were  taken 
to  determine  the  substrate  orientation  dependence  of  the 
c-axis  orientation.  The  {211}  orientation  was  determined 
to  suppress  all  but  one  c-axis  orientation  by  the  presence 
of  four  asymmetric  additional  spots  in  the  Laue  photograph. 
Therefore  the  {  211 }  orientation  of  the  GaAs  substrate 
best  suppresses  antiphase  domain  boundaries  associated 
with  material  regions  having  different  c-axes. 


The  second  reason  the  (211}  substrate  orientation 
is  important  for  this  work  is  related  to  the  suppression 
of  net  interfacial  charge.  Arguments  presented  by  W.A. 
Harrison  and  co-workers  [7]  show  that  certain  orientations 
of  polar/nonpolar  or  binary  on  elemental  semiconductor 
interfaces  will  tend  to  minimize  interfacial  charge. 

These  arguments  can  be  extended  to  ternary  on  binary 
interfaces.  The  development  of  these  concepts  is  somewhat 
involved,  so  here  we  present  a  simplified  argument  utilizing 
Harrison's  transmutation  alchemy.  Figure  7  shows  the  inter¬ 
face  between  a  hypothetical  ternary  compound  InGaP2  and  the 
binary  compound  GaAs.  The  orientation  of  the  interface 
is  shown  parallel  to  one  of  the  <111>  bond  directions.  Since 
each  plane  of  atoms  parallel  to  the  interface  in  this 
case  contains  equal  numbers  of  anions  and  cations,  these 
planes  are  on  average  charge  neutral.  Slight  misorienta- 
tions  of  this  interface  with  respect  to  the  <11 1>  bond 
directions  will  not  on  average  introduce  significant 
interfacial  charge  since  we  are  considering  III-III-V2  on 
III-V  system.  Harrison  points  out  that  the  tetrahedral 
bond  configuration  assures  that  each  of  the  bonds  connecting 
each  atom  to  its  four  nearest  neighbors  contains  exactly 
two  electrons.  The  electron  distribution  along  each  bond 
depends  on  which  species  are  being  bonded,  but  the  overall 
bond  charge  remains  the  same.  Using  this  argument  one 
can  determine  the  overall  interface  charge  by  counting 


Transfer  1  proton: 
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Figure  7 


Charge  neutral  interface 
construction. 


each  column  V  atom  as  having  +1  charge  from  the  extra  proton 
relative  to  column  IV,  and  each  column  III  atom  as  having 
-1  charge.  One  can  see  that  according  to  this  approach,  in 
this  III-III-V2  on  Ill-v  system,  small  misalignments  of  the 
interface  with  respect  to  the  <111>  bond  directions  will 
not  produce  significant  interfacial  charge. 

The  situation  is  very  different  for  mixed  column 
materials  systems  such  as  the  II-IV-V2  ZnSnP2  on  II I -V  GaAs 
he ter o junction.  To  examine  this  case,  imagine  the  hypo¬ 
thetical  transfer  of  a  proton  from  each  of  the  Ga  atoms  to 
each  of  the  In  atoms  as  in  Figure  7.  This  process  will 
transmute  Ga  into  Zn  and  In  into  Sn  thereby  producing 
ZnSnP2  on  GaAs.  Assuming  that  the  electron 
configuration  during  this  proton  transfer  remains 
unchanged,  the  only  charge  transfer  of  concern  is 
that  due  to  the  proton  movement.  This  ignores  the  part 
of  the  electron  distribution  that  would  tend  to  follow 
the  proton  motion.  Therefore  the  real  positive  charge 
transfer  for  each  proton  moved  is  less  than  one.  For 
simplicity  we  ignore  these  complications.  Since  the 
interface  is  in  this  case  parallel  to  the  <111>  bond 
directions,  the  positive  charge  moves  parallel  to  the  in¬ 
terface  and  so  the  interfacial  plane  remains  on  average 
charge  neutral.  If  we  misalign  the  interface  from  the 
<111>  bond  directions,  however,  there  will  be  a  component 
of  the  positive  charge  motion  either  toward  or  away  from 
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the  interface.  This  will  result  in  the  accumulation 
of  positive  or  negative  charge,  respectively,  at  the  in¬ 
terface.  The  charge  motion  in  the  other  parallel  planes 
of  the  top  layer  will  on  average  cancel  in  pairs. 

A  mathematical  way  of  expressing  the  condition  for 
a  charge  neutral  interface  can  be  stated  as  a  relation 
between  the  Miller  indices  (hkl)  of  the  interfacial  plane 
and  the  <111>  bond  direction.  Let  [hkl]  indicate  the 
direction  perpendicular  to  the  interfacial  plane.  For 
the  interfacial  plane  to  be  parallel  to  the  <111>  bond 
directions,  the  [hkl]  direction  must  be  perpendicular  to 
that  direction.  This  gives  the  condition, 

[hkl]  •  <111>  *  ±h  *  k  ±  1  »  0  .  (2.1) 

The  lowest  index  orientation  that  can  satisfy  this 
condition  is  the  {110}  orientation.  This  is  a  happy 
coincidence  for  some  applications  since  the  {110}  planes 
are  the  natural  cleavage  planes  for  III-V  compounds. 
Unfortunately  in  our  case,  growth  of  ZnSnP2  on  {110}GaAs 
is  not  the  best  for  the  elimination  of  antiphase  domain 
boundaries.  The  next  lowest  index  that  satisfies  the 
charge  neutral  condition,  the  {211}  orientation,  is  also 
that  which  best  suppresses  the  formation  of  these  antiphase 
domain  boundaries.  These  are  the  two  important  aspects  of 
the  {211}  substrate  orientation  for  our  devices. 
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Interfacial  charge  is  a  problem  in  hetero junctions 
because  among  other  effects,  it  can  greatly  alter  the 
experimental  determination  of  the  band  lineup.  Although 
interfacial  dipoles  would  be  necessary  to  actually  change 
the  band  offsets,  interface  charge  will  affect  the  experi¬ 
mental  values  of  the  band  offsets  derived  from  C-V  data. 
These  complications  will  be  discussed  further  in 
Sections  2.1.4  and  2.3. 

2.1.2  Absorption  Experiments 

Both  absorption  and  emission  experiments  were  performed 
at  room  temperature  on  these  ZnSnP2/GaAs  devices  to  test  our 
model  of  the  below  bandgap  mechanism  for  absorption  and 
emission.  First  we  shall  consider  the  results  of  the 
absorption  experiments.  The  apparatus  used  for  these 
experiments  is  shown  in  Figure  8 .  The  source  of  mono¬ 
chromatic  radiation  incident  on  the  hetero junction  was 
a  halogen  vapor  lamp  focused  by  a  toroidal  mirror  on  the 
entrance  slit  of  a  Spex  Minimate  monochromator.  The 
entrance  and  exit  slits  of  the  monochromator  were  set 
at  1.25mm  to  maximize  resolving  power  while  still  main¬ 
taining  an  acceptable  signal  level.  This  slit  width 
results  in  a  50  A  spectral  bandpass.  The  exit  radiation 
was  chopped  at  700  Hz  and  was  then  incident  on  the  cleaved 
face  of  the  sample.  This  cleaved  plane  is  perpendicular 
to  the  hetero-interface.  For  the  photovoltaic 
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Figure  8  Absorption  experiment  apparatus 


experiments,  no  voltage  bias  was  applied  to  the  sample. 

The  ac  response  voltage  across  the  sample  was  measured 
directly  by  an  EG&G  12 4A  lock-in  amplifier.  The  phase 
was  adjusted  to  maximize  the  signal  and  the  time  constant 
was  set  to  as  low  a  value  as  possible  while  avoiding 
excess  noise.  The  wavelength  was  scanned  automatically 
from  0.7  pm  to  1.1pm  by  the  wavelength  drive  motor.  The 
controlling  voltage  to  the  wavelength  drive  served  as 
the  x-axis  of  an  HP7035B  x-y  recorder.  The  chopper  motor 
served  as  the  external  frequency  reference  of  the  lock-in 
amplifier.  The  lock-in  output  was  used  as  the  y-axis 
of  the  x-y  recorder.  The  wavelength  was  scanned  in  the 
slow  scan  mode  which  took  eight  minutes  to  scan  /u  =  0.4pm. 

The  photovoltaic  response  as  a  function  of  photon 
energy  is  shown  in  Figure  9.  The  scale  of  the  x-axis 
is  actually  linear  in  wavelength  and  not  energy,  but  is 
shown  in  energy  units  for  ease  of  interpretation.  Two 
broad  peaks  are  observed  in  the  room  temperature  photo¬ 
voltaic  response.  The  larger  peak,  having  a  half -power 
point  at  1.42eV,  appears  to  be  the  band-to-band  absorption 
peak  for  GaAs.  Since  the  ZnSr^  is  much  more  heavily  doped 
than  the  GaAs,  most  of  the  depletion  region  occurs  in 
the  GaAs.  In  the  absorption  mode,  the  signal  comes  only 
from  those  carriers  created  where  there  is  a  field  to 


separate  them.  Therefore,  the  detectable  absorption  occurs 
only  in  the  depletion  region  and  so  it  is  reasonable  that 
no  ZnSnP2  absorption  peak  is  observed.  The  smaller  peak 
has  its  maximum  at  1.36eV.  The  amplitude  of  this  peak 
depends  on  the  orientation  of  the  he tero- interface 
to  the  incident  radiation  as  can  be  seen  in  Figure  10. 

The  below  bandgap  peak  almost  disappears  when  the  radiation 
is  incident  normal  to  the  epitaxial  layer  rather  than  to  the 
cleaved  face.  The  orientation  of  the  normal  to  the  cleaved 
face  with  respect  to  the  incoming  light  varied  from  0°  to 
45°  for  maximum  signal.  These  results  may  indicate  that 
the  below  gap  absorption  depends  on  the  polarization  of  the 
incident  light  with  respect  to  the  interface.  The  largest 
below  gap  response  occurs  when  the  polarization  is  in  a 
plane  perpendicular  to  the  interface,  or  when  the  light  is 
incident  parallel  to  the  interface. 

If  we  assume  that  this  below  bandgap  absorption  peak 
is  due  to  the  mechanism  that  has  been  described,  one  can 
extract  an  estimate  of  the  interfacial  energy  gap  Ej  from 
the  absorption  data.  Looking  back  at  Figure  3,  it  is 
evident  that  no  significant  absorption  will  take  place 
at  energies  below  the  interfacial  energy  gap  Ej  since  this 
would  involve  highly  improbable  tunneling.  Therefore 
to  extract  E^  from  the  below  gap  peak  it  is  not  appropriate 
to  use  either  the  maximum  or  half -power  point  energies. 
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Figure  10  Orientation  dependence  of  illumina¬ 
tion  on  below  bandgap  photoresponse. 


To  estimate  E  ^  roughly  we  can  extrapolate  the  steepest 
part  of  the  below,  gap  peak  to  the  zero  absorption  line. 

This  results  in  an  approximation  of  Ej  =  1.3 lev.  It  appears 
that  Ej  may  be  even  smaller  than  this  since  absorption  is 
observed  below  this  energy. 

We  also  performed  absorption  measurements  with  the 
sample  in  reverse  bias.  A  dc  voltage  was  applied  to  the 
sample  in  series  with  a  load  resistor.  The  signal  was 
taken  directly  off  the  series  resistor  and  fed  into  the 
lock-in  amplifier.  These  results  are  shown  in  Figure  11 
as  a  family  of  curves  for  biases  from  0  to  -1.6  V  separ¬ 
ated  by  equal  voltage  increments  (0.4V).  Both  the  high 
energy  absorption  and  the  below  gap  peak  increase  with 
increasing  reverse  bias.  The  increase  in  the  high  energy 
response  may  be  due  to  an  Increase  in  the  tunneling  of 
electrons  through  the  conduction  band  notch  into  the  GaAs 
where  the  electrons  were  excited  from  the  valence  band 
of  ZnSnP2.  This  explanation  is  reasonable  since  tunneling 
increases  with  field  and  so  reverse  bias.  The  below  bandgap 
absorption  should  also  increase  with  increased  field  in 
the  interfacial  region.  There  is  also  a  shifting  of  the 
peak  of  the  below  gap  response  with  reverse  bias.  This 
may  be  due  to  the  change  in  the  quasi-Fermi  levels  with 
reverse  bias.  As  reverse  bias  is  increased,  minority 
carriers  are  extracted  from  the  interfacial  region.  This 
results  in  fewer  electrons  or  filled  states  in  the 


conduction  band  well  of  ZnSnP2  and  fewer  holes  in  the 
valence  band  well  of  GaAs.  This  allows  absorption  at 
lower  energies  and  so  shifts  the  below  bandgap  peak  to 
longer  wavelengths. 

2.1.3  Emission  Experiments 

Next  we  conducted  room  temperature  emission  experiments 
on  these  ZnSnP2/GaAs  he tero junctions.  A  block  diagram  of 
the  apparatus  used  in  these  experiments  is  shown  in  Figure 
12.  This  apparatus  was  the  product  of  an  evolution  of 
attempts  to  measure  the  emission  spectra  of  these  devices. 
Short  range  emission  measurements  without  the  monochromator 
were  made  with  a  PbS  photoconductive  cell  held  at  liquid 
N2  temperature.  The  PbS  cell  was  selected  since  it  has 
detection  capability  at  wavelengths  greater  than  4ym. 

A  GaAs  filter  placed  in  front  of  the  PbS  cell  indicated 
that  there  was  significant  emission  below  the  absorption 
edge  of  GaAs.  The  PbS  cell  appeared  to  be  too  insensitive 
to  pick  up  the  emission  when  filtered  through  the  monochro¬ 
mator.  A  Si  p-i-n  photodiode  also  proved  to  be  too  insensi 
tive.  The  emission  from  the  hetero junction  with  constant 

forward  current  bias  of  20  to  100mA  or  current  density 

2 

of  10  to  50A/cm  could  be  seen  with  an  infrared  (IR)  viewer 
The  IR  viewer  was  used  to  line  up  the  hetero junction 
emission  onto  the  entrance  slit  of  the  monochromator. 


The  last  detector  used  was  a  Si  avalanche  photo  diode 
(APD).  Despite  the  fact  that  the  Si  absorption  cuts  off  at  a 
fairly  short  wavelength,  around  1  um  ,  the  increased  sensitivity 
of  this  detector  made  it  the  best  available  for  our 
experiments.  The  apparatus  for  the  emission  experiments 
was  then  set  up  using  the  Si  APD  as  shown  in  Figure  12. 

The  hetero junction  was  positioned  in  place  of  the  halogen 
vapor  lamp  using  a  two  degree  of  freedom  micromanipulator. 

The  hetero junction  emission  was  focused  with  a  toroidal 
mirror  onto  the  entrance  slit  of  the  monochromator.  There 
was  some  orientation  dependence  of  the  emission  with  the 
largest  signal  detected  when  the  normal  to  the  epitaxial 
layer  pointed  directly  into  the  toroidal  mirror.  The 
exit  emission  from  the  monochromator  was  chopped  at  700  Hz 
and  then  focused  onto  the  Si  APD.  The  Si  APD  was  biased 
appropriately  just  into  avalanche.  The  ac  signal  across 
a  47kn  series  resistor  housed  in  a  shielded  box  was  fed 
into  the  lock-in  amplifier.  The  lock-in  amplifier  was 
operated  in  the  differential  mode  with  careful  grounding 
of  the  detection  circuitry  to  ensure  maximum  signal  to 
noise  ratio.  The  calibration  of  the  monochromator  and 
alignment  of  the  Si  APD  was  checked  using  a  He-Ne  Laser. 

A  set  of  emission  curves  for  one  of  the  hetero junction 
emitters  is  shown  in  Figure  13.  The  sample  was  biased  at 
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Figure  13  ZnSnP2/GaAs  emission  spectra 
at  four  constant  current 
densities . 
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constant  current  densities  of  10,  20,  30,  and  40A/cm  . 

Two  peaks  are  observed  in  the  emission.  Consider  the 

2 

emission  peaks  at  the  lowest  current  bias  of  lOA/cm  . 

Here  the  slightly  larger  peak  has  its  maximum  at  1.42ev. 
This  peak  appears  to  be  the  band-to-band  emission  peak 
for  GaAs.  We  see  no  analogous  emission  peak  for  ZnSnP^. 
This  can  be  explained  in  a  similar  fashion  as  for  the 
absence  of  a  ZnSnP2  peak  in  the  absorption.  Since  the 
ZnSnP^  is  much  more  heavily  doped  than  the  GaAs,  many  more 
holes  are  injected  into  the  GaAs  well  than  electrons  into 
the  ZnSnP2  well.  Therefore,  there  are  many  more  minority 
carriers  available  in  the  GaAs  which  are  necessary  for 
the  band-to-band  radiative  recombination. 

The  smaller  emission  peak  has  its  maximum  at  1.36ev. 

If  we  again  assume  our  below  gap  emission  mechanism  is 
correct,  an  estimate  of  Ej  can  be  extracted  from  this  peak. 
As  in  the  case  for  absorption,  one  would  not  expect  to 
see  significant  emission  below  if  the  forward  bias 
diagram  is  correct.  We  can,  therefore,  estimate  Ej  by 
extrapolating  the  steepest  part  of  the  below  gap  emission 
peak  to  the  zero  emission  line.  This  results  in  an  esti¬ 
mate  of  Ej  =  1.3 lev  which  again  is  an  upper  limit.  This 
is  the  same  value  of  E^  as  that  obtained  from  the  absorp¬ 


tion  data. 
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Some  other  points  of  interest  concerning  these  emission 
data  should  be  mentioned.  One  notices  immediately  that 
the  two  emission  peaks  are  much  closer  in  height  than 
those  of  the  absorption  data.  In  the  case  of  absorption, 
the  active  region  for  the  band-to-band  GaAs  absorption 
is  the  depletion  region.  Using  the  depletion  approximation 
the  depletion  width  was  estimated  at  500 A.  The  active 
region  for  the  below  gap  absorption  is  much  smaller  than 
this  because  it  is  the  width  of  the  region  across  which 
significant  tunneling  can  occur.  This  is  likely  to  be 
a  few  lattice  constants  in  width  or  a  few  tens  of  A.  One 
would  therefore  expect  the  below  gap  absorption  peak  to 
be  approximately  an  order  of  magnitude  smaller  than  the 
band-to-band  absorption  peak.  In  the  case  of  emission, 
the  below  gap  and  band-to-band  mechanisms  are  both  con¬ 
trolled  by  the  supply  of  holes  injected  into  the  GaAs 
well.  These  holes  are  involved  in  the  radiative  re¬ 
combination  of  both  mechanisms.  Therefore,  it  is  reason¬ 
able  that  the  emission  peaks  are  closer  in  height  than 
those  of  the  absorption  spectra. 

Another  point  of  interest  is  the  shift  of  the  two 
emission  peaks  to  longr.r  wavelengths  at  increased  current 
bias.  The  most  likely  cause  of  this  shift  is  heating 
of  the  sample  since  the  bandgaps  of  both  semiconductors 
decrease  with  increasing  temperature.  The  shrinking  of  the 


bandgap  corresponds  to  an  increase  of  the  lattice  constant. 
Double-crystal  x-ray  diffraction  measurements  indicate  no 
difference  in  the  thermal  expansion  between  ZnSnP2  and 
GaAs  from  room  temperature  to  about  150 °C.  Since  the 
lattice  constants  on  both  sides  of  the  interface  expand 
together  at  the  same  rate,  we  expect  to  see  the  two 
emission  peaks  move  at  about  the  same  rate  with  respect 
to  current,  as  seen  in  Figure  13.  One  can  estimate  the 
temperature  changes  necessary  for  the  apparent  energy 
gap  variations  using  the  approximate  relation. 


Eg ( T )  =  E  (0) 


T+S 


(2.2) 


For  GaAs,  the  empirically  determined  constants  are 

a  =  5.405  x  10-4,  0  =  204,  and  Eg(0)  =  1.519eV  [8J.  For 
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a  constant  current  density  20  A/cm  ,  the  energy  gap  of 
GaAs  shifted  from  the  300K  value  of  1.42eV  to  1.41eV. 
Inserting  AEg=0.01eV  into  equation  (2.2)  and  solving  the 
quadratic  yields  AT  =  27.4K.  So  for  all  values  of  current 
density  used  in  these  experiments,  the  temperature  of 
the  device  should  remain  under  150 °C.  Therefore,  the 
behavior  of  the  emission  peaks  is  consistent  with  that 
expected  from  device  heating. 

Now  compare  the  absorption  and  emission  results 


directly  by  examining  Figure  14.  Here  the  higher  solid 
curve  is  the  photovoltaic  response  whereas  the  lower 
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Figure  14  Overlay  of  photovoltaic  response 
and  low  bias  emission  for 
ZnSnP2/GaAs. 
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dashed  curve  is  the  emission  spectrum  at  the  lowest  current 
2 

bias  of  10  A/cm  .  The  higher  energy  emission  peak  has  its 
maximum  at  the  half -power  point  energy  of  the  larger  absorp¬ 
tion  peak.  One  expects  this  type  of  relation  between 
absorption  and  emission  for  a  band-to-band  transition. 

The  energies  of  these  peaks  are  appropriate  for  GaAs. 

The  below  bandgap  peaks  in  the  emission  and  absorption 
spectra  are  coincident  in  energy.  This  observation  may 
help  to  address  an  obvious  concern  in  analyzing  these 
data.  There  is  the  possibility  that  the  below  gap  peaks 
are  not  due  to  the  mechanisms  described  but  rather  to 
transitions  involving  tunneling  to  and  from  impurity  levels 
in  the  gap.  The  energy  coincidence  of  the  two  below  gap 
peaks  might  be  taken  as  evidence  against  an  impurity  inter¬ 
pretation.  When  optical  transitions  of  electrons  involve 
impurity  levels,  in  principle  one  expects  a  shift  between 
the  absorption  and  emission  peaks,  where  the  emission  peaks 
at  longer  wavelengths.  This  shift  is  due  to  the  Franck- 
Condon  principle  [9,10]  which  states  that  electron 
transitions  are  so  fast  that  the  lattice  can  be  viewed 
as  static  during  the  transition.  These  issues  will  be 
addressed  in  Section  2.4. 

2.1.4  Band  Lineup 

The  analysis  of  the  emission  and  absorption  data 
in  Sections  2.1.2  and  2.1.3  has  been  based  on  the  assumption 


that  the  band  lineup  of  p-ZnSnP2/n-GaAs  is  staggered.  To 
show  that  the  lineup  is  indeed  staggered  and  so  extract 
a  theoretical  value  of  ,  we  need  to  measure  several 
parameters  of  the  materials  themselves  as  well  as  of  the 
hetero junction  as  a  whole. 

To  construct  the  band  lineup,  we  first  determined 
the  built-in  potential  V  .  The  built-in  potential  is 

the  difference  between  the  chemical  potentials  of 
the  two  materials  or  the  potential  that  arises  across 
the  "hetero junction  from  charge  that  transfers  to  equilibrate 
the  electrochemical  potential  throughout  the  device.  We 
used  the  C-V  intercept  method  to  determine  V  . 

For  an  abrupt  hetero junction  of  area  A  and  doping 
Na  and  permittivity  ep  for  the  p-type  material  and 
Nd  and  en  for  the  n-type  material,  one  can  derive  the 
relation. 


c-2  ,  2(SWa  +  W(VV> 

qenepA2NdNa 

This  equation  reveals  a  linear  relationship  between  C 


(2.3) 

-2 


and  the  applied  potential  V  [11].  The  built-in  potential 

-2  -2 

VQ  is  obtained  by  extrapolating  the  C  -V  plot  to  C  =0 
where  the  voltage  intercept  is  V  =  V  .  This  expression 
neglects  the  effects  of  interface  charge  Interface  charge 
effects  will  be  discussed  later  in  this  section,  and  in 
Section  2.3. 
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Capacitance  versus  voltage  measurements  were  performed 

using  a  Boonton  7  2BD  capacitance  meter.  The  Boonton  meter 

operates  at  1  MHz.  The  capacitance  of  the  leads  and  sample 

holder  were  nulled  out  using  a  nulling  circuit  in  the 

Boonton  meter.  The  sample  was  reverse  biased  out  to 

-2  volts,  beyond  which  the  sample  went  into  breakdown. 

-2 

A  plot  of  C  versus  V  for  these  samples  is  shown  m 

Figure  15.  The  slope  and  the  intercept  were  determined 

using  linear  regression  analysis.  The  built-in  potential 

was  determined  as  =  1.75  V. 

o 

Before  proceeding  to  construct  the  band  lineup,  let 
us  compare  this  measured  value  of  VQ  with  a  calculated 
value.  As  mentioned  earlier,  the  built-in  potential  is 
the  difference  between  the  chemical  potentials  un  and 
Up  in  the  two  materials, 

«Vo  *  “n  *  up  •  <2'4) 

Adding  and  subtracting  the  band  edge  energies  yields, 

«Vo  =  -  Ec  >  +  Ec  ~Ev  -  VEv  (2'5) 

n  n  p  *  p 

Since  both  semiconductors  are  degenerately  doped,  to  de¬ 
termine  the  terms  in  parentheses  we  must  use  the  following 
expressions  for  nondilute  solutions, 


_  2 

Figure  15  2nSnP2/GaAs  C  -V  plot  shows 

V  a  1.75  V  at  1  MHz. 
o 
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and 
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'P 


/%-M 

*  \ F*  \~V— A 


(2.7) 


where 


W 


f  n^dn  _ 

1  +exp(n-  Hr ) 
o  1 


(2.8) 


and 


/2irm  *kT\  3/2 
N<:n  *  ^  \  /  • 


(2.9) 


Inverting  equations  (2.6)  and  ( 2 . 7 ),  equation  (2.5) 
becomes , 


’vo  =  Ec„  -  Ev„  +  kT 
n  p 


f4*1(5!1  1  +  f»"1(5E  1 

*  n<3n  *  V 


(2.10) 


Now  if  one  assumes  for  this  calculation  that  the 
lineup  is  staggered,  then  a  glance  at  the  flat  band  diagram 
of  Figure  2  shows  that, 


and 


E  -  E  =  AE  +  E  , 
c_  v  v  g  ' 
n  p  3n 


AE„  =  AE_  +  A  E_. 
vc  g 


(2.11) 


Before  making  a  determination  of  aEv,  we  should  mention 
some  of  the  controversies  associated  with  determinations 
of  band  offsets.  One  of  the  oldest  methods  for  determining 
the  conduction  band  offset  is  through  the  use  of  the  electron 
affinity  rule  [12].  This  states  that  aE  is  equal  to  the 


„V 


*3 


a 


.V 

i 


'-3 


difference  in  the  electron  affinities  of  the  two 


semiconductors , 

AEC  =  X1_X2  =  AX*  (2.12) 

This  method  is  plagued  with  fundamental  and  practical 
problems.  On  the  practical  side,  one  is  attempting  to 
determine  a  relatively  small  number,  aE  ,  on  the  order  of 
a  few  tenths  of  an  eV  from  two  much  larger  numbers,  x-^  and 
x2»  which  are  typically  a  few  eV  in  size.  If  there  is 
any  significant  error  in  the  values  of  x^  and  x2#  AEC  will 
have  an  error  on  the  order  of  its  own  magnitude.  Also, 
there  is  evidence  that  the  electron  affinity  values  (x ) 
may  depend  on  the  crystallographic  plane  for  which  they 
were  measured.  The  x  for  a  particular  semiconductor  for 
the  proper  orientation  may  not  be  available. 

Problems  exist  with  this  "rule"  on  a  more  fundamental 
level.  The  electron  affinity  is  essentially  a  property 
encompassing  aspects  of  both  the  bulk  of  the  semiconductor 
and  its  surface.  When  an  electron  passes  from  one  semi¬ 
conductor  to  the  other  in  a  hetero junction,  however,  the 
surface  effects  are  eliminated  and  replaced  by  interface 
effects.  So  we  chose  to  approach  the  band  offset  determina¬ 
tion  by  considering  the  bulk  rather  than  the  surface  proper¬ 
ties  of  the  two  semiconductors.  We  approximated  the  valence 
band  offset  as  the  difference  in  the  ionization  energies  of 
the  semiconductors, 


A  E. 


Al. 


(2.13) 
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This  assumes  that  no  interfacial  dipole  is  formed  when 
the  two  materials  are  joined.  The  ionization  energy 
is  defined  as  the  energy  difference  between  the  top  of 
the  valence  band  and  the  vacuum  state.  The  measured 
ionization  energy  is  related  to  surface  properties  and 
is  usually  only  available  for  natural  cleavage  planes 
for  each  semiconductor.  To  approach  the  ionization  energy 
as  a  bulk  property,  we  derived  the  ionization  energies 
of  ZnSnP^  and  GaAs  using  a  dielectric  two-band  model 
[13,14,15].  The  results  of  that  calculation  are' 


IZnSnP2  =  5,91eV 


and 


LGaAs 


*  5.68eV. 


This  gives  AEv=0.23eV.  To  complete  the  calculation 
of  V  ,  we  need  to  insert  values  for  the  carrier  concentra¬ 
tions  and  the  appropriate  effective  densities  of  states. 

* 

Inserting  the  electron  effective  mass  for  GaAs,  m  = 0.0665m 

C  \ 

[16],  into  equation  (2.9)  for  the  effective  density  of 

states  one  obtains  for  GaAs,  Nc  =  4.3  x  1017cm"3. 

An  approximation  for  the  hole  effective  mass  for 

ZnSnP2  is  necessary  since  no  experimentally  determined 

values  are  currently  available.  Using  a  theoretical  value 

* 

determined  by  Goryunova  and  co-workers  [17]  of  =  0.39mQ, 

the  effective  density  of  states  of  the  valence  band  of 

18  -  3 

ZnSnP2  is  calculated  as,  N  =  6.1  x  10  cm 

P 


i 
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Using  the  carrier  concentrations  determined  by  Hall 

18  -3  19  -3 

measurements  of  n_ .  =  10  cm  and  pn  =  5  xlO  cm  ,  and 

n  p 

the  Fermi  Dirac  integrals  from  tables,  the  built-in 
potential  is  calculated  as,  V  «  1.81  V.  This 
is  in  good  agreement  with  the  value  obtained  from  the 
C-v  data  (1.75V),  considering  the  approximations  made 
in  this  calculation. 

There  is  a  fundamental  problem  with  the  C-v  intercept 

method  of  determining  V  .  Equation  (2.3)  between  and 

V  neglects  the  effects  of  interfacial  charge.  An 

expression  for  C  which  is  valid  in  the  presence  of 

voltage  independent  interface  charge  has  been  derived  [18]. 

This  analysis  shows  that  interfacial  charge  reduces 

the  measured  value  of  VQ  regardless  of  the 

sign  of  the  charge.  All  of  our  C-V  measurements  for 

_  2 

ZnSn^/GaAs  devices  yielded  linear  C  -V  plots.  Earlier 

researchers  presumed  this  linearity  ensured  that  interfacial 

charge  was  insignificant.  In  fact,  interfacial  charge  will 

shift  the  intercept  V  =  VQ  to  smaller  values  but  will  not 

affect  the  slope  of  the  line.  Based  on  this  information, 

we  believe'  the  most  accurate  measured  value  of  V^  would  be 

the  largest  one  we  obtained.  The  measured  VQ  values  ranged 

from  1.33  to  1.75  V.  All  V  values  obtained  from  the  C-v 

o 

data  except  the  smallest  value  of  1.33  V,  confirm  the 
staggered  lineup  of  p-ZnSnl^/n-GaAs .  Despite  the  fact 
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that  ideally  we  are  dealing  with  a  charge  neutral  orienta¬ 
tion,  real  substrate  orientations  can  only  be  prepared 
accurately  to  ±0.5°.  This  small  misorientation  can  result 
in  significant  interfacial  charge  in  our  hetero junctions. 
Considering  this  we  shall  use  the  value  of  V  =  1.75  V  for 
our  band  lineup  calculation. 


We  now  calculate  the  band  lineup  of  p-ZnSnP2/n-GaAs  to 
show  that  it  is  indeed  staggered  and  to  provide  a  prediction 
for  the  interfacial  energy  gap  E_.  Looking  back  at  Figure  2 
(the  flat  band  diagram  for  this  hetero junction)  one  sees 
that  Ej  can  be  calculated  if  we  obtain  the  band  offsets  and 
the  separation  of  the  appropriate  band  edges  from  the 
chemical  potentials  in  the  two  materials. 

Inverting  equations  (2.6)  and  (2.7)  for  the  degenerate 


carrier  concentrations  we  obtain. 


From  the  flat  band  diagram  it  is  apparent  that, 


iE  =  qV  -  ( u  -  E  )  -  E  -  ( E  -  u  )  . 
v  M  o  Hn  c  g„  v  p 

n  n  p 

This  results  in  AEV  =  0.175  eV.  Also  we  find 

AE  =  AE  -AE  = 
c  v  g 


0.145  eV. 
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Now  the  interfacial  energy  gap  can  be  calculated  as, 

E_  =  a V  -  (u  -  E  )  -  AE  -  A E  -  (E  -  u  )  . 

I  H  o  VMn  c  c  v  v_ 

n  p 

This  results  in  E-  *  1.275  eV.  This  is  also  in  good 

xTh 

agreement  with  the  experimental  value  of  ET  *  1.31  eV. 

1EX 

As  mentioned  before,  the  true  experimental  value  of  Ej 
is  probably  even  smaller  than  that  indicated  by  the  extra¬ 
polation  of  the  data  and  so  the  theoretical  value  may 
be  more  accurate  than  it  appears.  Still,  it  is  necessary 
to  address  the  problems  in  our  interpretations  of  the 
data,  since  other  explanations  cannot  be  ruled  out.  These 
will  be  discussed  in  the  conclusions  section. 

2.1.5  Other  Results 

Other  experiments  were  performed  on  these  devices 
to  understand  various  properties  of  the  hetero-interface. 
These  efforts  were  primarily  aimed  at  measuring  device 
capacitance  as  a  function  of  frequency  and  illumination. 
Here  the  procedures  and  results  of  these  experiments  are 
discussed  along  with  suggestions  for  improvements  for 
future  investigations. 

There  are  several  nondestructive  electrical  techniques 
involving  the  measurement  of  capacitance  to  obtain  infor¬ 
mation  about  impurity  concentrations  and  trap  levels  in 
semiconductor  devices  [19].  For  instance,  a  small 


ac  voltage  applied  to  a  device  can  cause  the  charging 
and  discharging  of  trap  levels,  thereby  affecting  the 
device  capacitance.  One  popular  technique,  used  in  the 
analysis  of  an  oxide-semiconductor  interface  involves 
sweeping  the  Fermi  level  past  interface  state  energy 
levels  by  changing  the  applied  bias.  This  will  cause 
the  occupancy  of  the  trap  state  to  change  and  so  will 
contribute  to  the  capacitance  and  conductance  of  an  MOS 
device.  Also,  the  charging  and  discharging  of  interface 
states  may  be  rate  limited  because  of  their  relaxation 
times.  This  will  result  in  a  frequency  dependence  of 
the  capacitance  associated  with  interface  states.  We 
reasoned  that  these  ideas  and  techniques  might  also  be 
applicable  to  the  characterization  of  interface  states 
and  charging  mechanisms  in  a  he ter o junction.  These 
different  mechanisms  may  have  different  frequency 
dependences  which  we  hoped  to  decipher.  Since  the 
occupation  of  interface  states  as  well  as  the  number  of 
free  carriers  can  be  affected  by  illuminating  the  sample, 
we  also  conducted  photocapacitive  experiments  at 
different  frequencies. 

Before  setting  up  for  variable  frequency  capacitance 
experiments,  we  first  performed  simple  experiments  to  see 
the  general  effect  of  light  on  device  capacitance  at 
discrete  frequencies  using  available  capacitance  meters. 
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First  using  an  HP4261A  LCR  meter,  we  measured  the  change  in 
device  capacitance  at  1  kHz  as  a  function  of  the  incident 
wavelength  from  0.7  to  1.0  yin.  The  results  are  shown  in 
Figure  16.  These  data  were  obtained  by  shining  unchopped 
monochromatic  light  on  the  sample  and  reading  the  capaci¬ 
tance  directly  off  the  HP  meter.  Since  application  of  dc 
light  bias  is  very  similar  to  the  application  of  dc  voltage 
bias,  the  HP  meter  had  to  be  used  in  the  configuration 
appropriate  for  C-V  measurements.  This  entailed  the  in¬ 
sertion  of  a  1.1  yF  blocking  capacitor  between  the  high 
current  and  voltage  terminals.  The  HP  meter  cannot  pro¬ 
cess  a  dc  current  as  was  easily  verified  by  shorting  these 
terminals  together.  In  this  configuration,  the  dark  capa¬ 
citance  was  correctly  determined  but  no  change  in  C  was 
registered  when  the  sample  was  illuminated. 

The  wavelength  was  varied  in  increments  of  10  nm  using 
a  slit  width  of  1.25  mm  which  limits  the  resolution  to 
±2.5  nm.  At  1  kHz  the  capacitance  appears  to  have  a  spec¬ 
tral  variation  very  similar  to  that  of  the  absorption  data 
with  peaks  at  1.44  eV  and  1.36  eV.  To  extract  the  capaci¬ 
tance  of  the  device,  the  HP  meter  modelled  these  samples 
as  an  ideal  capacitance  in  parallel  with  an  ideal  resis¬ 
tance.  The  meter  determined  and  displayed  the  capacitance 
and  dissipation  factor.  For  a  parallel  RC  circuit,  the 
dissipation  factor  is  equal  to  D  =  l/(coRpCp).  From  D  at 


Wavelength  (Mm) 


Figure  16  Capacitance  as  function  of 

photon  energy  of  ZnSnP2/GaAs 
at  1  kHz . 


each  wavelength,  the  conductance  Gp  =  1/Rp  was  obtained  and 
is  plotted  in  Figure  17.  The  conductance  and  capacitance 
spectra  are  very  similar  except  that  the  conductance  peaks 
are  broader.  Capacitance  and  conductance  spectra  at  a 
single  frequency  such  as  these  could  not  provide  concrete 
evidence  to  support  or  refute  our  below  bandgap  model. 

We  needed  therefore  to  explore  these  spectra  at  many  fre¬ 
quencies  . 

The  HP  meter  also  has  a  120  Hz  frequency  setting. 

The  capacitance  of  these  devices  at  120  Hz  fluctuated 
wildly  and  the  dissipation  factor  was  off-scale.  The 
expression  for  the  dissipation  factor  reveals  the  problem. 
At  low  frequencies  the  dissipation  factor  becomes  too 
large  for  the  measurement  of  C  to  be  accurate.  To  resolve 
the  problem,  a  large  parallel  capacitance  was  added  across 
the  sample.  The  variation  of  C  with  X  at  120  Hz  was  deter¬ 
mined  in  this  manner  to  be  roughly  the  same  as  that  at 
1  kHz . 

A  Boonton  capacitance  meter  was  used  to  perform 
capacitance  measurements  at  1  MHz.  Very  little  change 
in  the  capacitance  with  illumination  was  observed.  This 
result  is  probably  due  all  or  in  part  to  the  measurement 
configuration  of  the  Boonton  meter  rather  than  the  sample 
itself.  The  Boonton  meter  is  housed  in  a  cabinet  with  the 
necessary  electronics  to  perform  C-V  measurements.  The 
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Figure  17  Photoconductance  at  1  kHz  for 

ZnSnP2/GaAs  taken  from  dissipa¬ 
tion  measurements. 
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device  capacitance  is  measured  by  the  Booton  with  a  two 
terminal  configuration.  Unless  some  modifications  are 
made,  this  configuration  probably  will  not  detect  changes 
in  the  capacitance  which  contribute  to  a  dc  current. 

To  measure  the  capacitance  of  these  devices  over  a 
wide  frequency  range  both  with  and  without  illumination, 
we  chose  to  implement  a  technique  using  a  lock-in  ampli¬ 
fier.  The  EG&G  lock-in  has  its  own  internal  reference 
oscillator  which  can  generate  sine  waves  of  frequencies 
from  0.2  Hz  to  210  kHz.  The  capacitance  measurement  con¬ 
figuration  is  shown  in  Figure  18.  The  internal  ac  refer- 
voltage  vT  was  applied  across  the  series  combination  of 
the  sample  and  a  variable  load  resistor.  The  lock-in 
amplifier  was  used  in  the  differential  mode  measuring  the 
difference  of  the  input  terminal  voltages,  vA  -  vB.  The 
reference  voltage  was  applied  across  the  B  input  while  the 
sample  voltage  served  as  the  A  input.  The  resulting  total 
input  to  the  lock-in  was  the  voltage  vR  across  the  load 
resistor.  By  recording  the  0°  or  in  phase  signal  voltage 
v0 ,  the  90°  signal  voltage  vg0,  and  the  total  reference 
voltage  vT,  the  parallel  capacitance  and  parallel  resis¬ 
tance  of  the  sample  were  determined  at  each  frequency. 

The  reference  voltage  was  measured  by  using  the  lock-in 
in  its  ac  voltmeter  mode. 
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The  derivation  of  the  equations  used  to  determine 
the  capacitance  and  resistance  of  the  sample  is  somewhat 
lengthy,  so  here  an  outline  of  the  derivation  is  given. 
To  determine  the  voltage  drop  across  the  load  resistor, 
one  can  use  voltage  division  to  obtain. 


v„  = 


*L  VT 


R  RL+[l/Rp+jWCp: 


(2.15) 


For  simplicity,  let  Rp  =  R  and  Cp  =  C.  Multiplying  numer¬ 
ator  and  denominator  by  the  appropriate  terms  gives. 


[RL/R  +Rl2/R2+<*2R2lC2+  ju»RLc]  vT 
R  (1+Rl/R)2  +w2rl2c2 


(2.16) 


O 

The  0  or  in  phase  voltage  is  then; 


(R^R  +  R^/R2  +  ‘-2Rl2C2  Jv^ 
°  (1  +  Rj^/R)2  +  u)2R^2C2 


The  90  voltage  is, 


(2.17) 


(uRlC)vt 

r90  =  2  2  2  2  * 

(i  +  rl/R)  +  <*>  rl 


(2.18 


From  equation  (2.17), 


(1  +  Rl/R)  vt(Rl/R)  -  vq(R/Rl+1) 

C'L  w2rl2  (vq  -  vT) 


(2.19 
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From  equation  (2.18), 

(u2!^2  vgo)  C2  -  (wRIyT)  c  +v90(l  +  Rj/R)2  -  0.  (2.20) 

Using  equations  (2.19)  and  (2.20)  to  eliminate  C  and 
determine  the  unknown  parallel  resistance  R  we  obtain, 


vT(vrr  -  vQ ) 
+  l  -  "2 


(2.21) 


+  <VT-V0) 


Equation  2.20  can  be  solved  to  yield, 


The  appropriate  sign  of  the  square  root  in  this  expression 
is  selected  to  give  the  most  reasonable  value  of  c  at 


each  frequency. 


The  reference  voltage  was  set  at  15  mV  rms.  The 
load  resistor  was  varied  from  a  few  hundred  ohms  to  about 
50  kn.  Lower  values  of  R^  yielded  flatter  and  probably 
more  accurate  frequency  spectra  for  the  capacitance. 
Regardless  of  the  size  of  R^,  there  was  a  tendency  for  the 
capacitance  of  the  device  to  rise  at  low  frequencies  above 
its  plateau  value  of  about  600  pF.  At  high  frequencies 
the  measured  capacitance  dropped.  The  measured  resistance 
tended  to  remain  steady  at  lower  frequencies  and  then 
rapidly  drop  at  high  frequencies. 


It  appears  that  most  of  the  observed  capacitance 
and  resistance  changes  of  these  devices  were  due  to  in¬ 
accuracies  in  the  measurement  procedure.  This  conclusion 
was  reached  primarily  by  measuring  the  capacitance  and 
resistance  of  a  known  RC  circuit  as  a  function  of  fre¬ 
quency.  The  values  for  the  RC  circuit  were  chosen  at 
the  1  kHz  RC  values  of  the  sample  and  so  the  circuit  was 
called  the  mock  sample.  The  mock  sample  showed  much  of 
the  same  trends  as  a  function  of  frequency  as  the  real 
device.  One  possible  source  of  error  which  might  lead 
to  these  results  is  error  in  the  phase- shifting  circuitry 
of  the  lock-in.  This  error  is  indicated  in  the  lock-in 
manual  as  being  worst  at  very  low  and  very  high  frequencies. 
Without  accounting  for  this  error,  the  capacitance  and 
resistance  could  at  best  be  measured  accurately  only 
between  50  Hz  and  5  kHz. 

The  change  in  capacitance  and  resistance  of  these 
devices  upon  illumination  was  determined  using  the  lock-in 
from  20  Hz  to  100  kHz.  At  each  frequency,  the  capacitance 
changes  of  the  sample  were  determined  when  illuminated 
with  the  two  peak  absorption  wavelengths  A=0.86ym  and 
A  =0.91 um.  Since  the  lock-in  amplifier  can  filter  out 
dc  current,  a  blocking  capacitor  is  not  needed.  According 
to  the  lock-in  measurements,  the  capacitance  of  the  device 


illuminated  with  A  =  0.86ym  was  3%  higher  than  the  dark 
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capacitance  over  the  whole  frequency  range.  With  this 
measurement  technique,  the  accuracy  of  the  absolute  values 
of  the  capacitance  at  each  frequency  is  probably  lower 
than  that  of  the  change  in  capacitance  with  illumination. 
Still  there  must  be  other  problems  with  the  measurement 
technique  since,  when  these  samples  were  illuminated  under 
the  same  conditions  and  measured  with  the  HP  meter,  a 
10%  increase  in  capacitance  was  observed. 

Two  improvements  in  the  capacitance  measurement  using 
a  lock-in  amplifier  can  be  made.  Firstly,  the  ac  voltage 
across  the  sample  should  be  kept  at  the  same  value  for 
the  whole  frequency  range.  This  would  entail  adjusting 
the  series  load  resistor  or  reference  voltage  at  each 
frequency.  Secondly,  since  error  in  the  phase  determina¬ 
tion  of  even  0.1°  can  cause  great  inaccuracies  in  the 
measured  capacitance,  some  technique  for  assuring  proper 
measurement  of  the  phase  is  needed.  A  method  which  may 
be  applicable  here  is  one  which  was  used  to  measure  the 
forward-biased  capacitance  of  a  Schottky  barrier  diode 
[  19] .  The  problem  was  basically  one  of  measuring  a  small 
capacitance  in  parallel  with  a  large  conductance.  At 
very  low  frequencies,  the  in  phase  portion  of  the  current 
swamped  the  out  of  phase  current.  The  solution  proposed 
was  to  insert  a  small  resistor  in  parallel  with  the  sample 


If  considered  a  pure  resistance,  this  resistor  should 
not  affect  the  measured  value  of  the  device  capacitance 
at  any  frequency.  The  phase  was  adjusted  until  the 
capacitance  measured  with  and  without  the  test  resistor 
was  the  same.  In  our  case,  this  type  of  approach  may 
be  of  some  use,  but  proper  adjustment  of  the  ac  voltage 
across  the  sample  is  a  more  critical  issue. 


2.2  InjL_xGaxP/GaAs 

2.2.1  Material  Properties  and  Device  Construction 

As  mentioned  in  the  introduction,  it  has  been  pre¬ 
dicted  that  the  n (down) /p (up)  configuration  of  the  stag¬ 
gered  lineup  he tero junction  would  be  the  best  below  bandgap 
emitter  [2] .  To  achieve  this  configuration  in  the 
ZnSnP2/GaAs  system,  we  would  need  to  dope  the  ZnSnP2 
n- type .  Since  we  were  unable  to  do  this,  we  looked  to  a 
different  materials  system,  p-In^_xGaxP/n-GaAs .  This 
materials  system  has  more  flexibility  since  both  In^_xGaxP 
and  GaAs  can  be  made  n-type  or  p-type.  The  flat  band 
lineup  of  p-Ini_xGaxP/n-GaAs  is  shown  in  Figure  19.  The 
Ini-xGaxP  layers  were  grown  on  {lll}GaAs  substrates  using 
LPE  by  S. Julie  Hsieh  [3] .  Ini_xGaxP  is  a  III-V  ternary 
alloy  which  can  be  lattice-matched  to  GaAs.  To  determine 
the  band  lineup,  the  necessary  band  parameters  were  ob¬ 
tained  in  the  following  manner.  At  each  alloy  composition 
of  Ini„xGaxP,  the  energy  gap  was  determined  by  its  depen¬ 
dence  on  the  lattice  constant  which  was  obtained  by  double¬ 
crystal  x-ray  diffraction  scans.  The  lattice-matched  alloy 
In0. 49Ga0. 51p  was  determined  in  this  way  to  have  an 
energy  gap  of  1.90  eV.  The  energy  gap  of  the  GaAs  was 
measured  by  photoconductivity  on  high-purity  epitaxial 
samples  to  be  1.39  eV.  The  separation  of  the  chemical 


potentials  in  the  two  materials  from  the  closest  band  edge 
was  determined  by  Hall  measurements.  The  electron  concen¬ 
tration  of  the  GaAs  substrate  was  found  to  be 
n  *  1.4  x  10^®cm“3.  Hall  measurements  were  performed  on 
p-In,  Ga  P  layers  grown  on  SI  GaAs  under  the  same  condi- 
tions  as  those  grown  on  n-GaAs.  These  measurements  yielded  a 
hole  concentration  in  the  In1->xGaxP  of  p  =  2 . 5  x  1017cm”3 .  The 
built-in  potential  VQ  was  determined  by  C-V  measurements 
at  1  MHz. 

The  experimentally  determined  band  lineup  for  these 
{111}  p-In,  Ga  P/n-GaAs  hetero junctions  is  staggered,  but 
a  dielectric  two-band  calculation  predicted  the  lineup  to 
be  straddling.  Our  experimental  band  offsets  were 
AEV  *  0.08  eV  and  AEC  =  0.59  eV,  which  yields  Ex  =  1.31  eV. 
The  dielectric  model  predicts  AEV  =  -0.23  eV  and 
AEC  =  0.28  eV.  Note  that  the  {111}  interface  is  not,  allowing 
for  differences  in  electronegativity,  a  charge  neutral  inter¬ 
face.  We  conjecture  that  the  {111}  band  lineup  may  be 
shifted  from  the  predicted  straddling  lineup  to  a  staggered 
one  due  to  dipoles  associated  with  interface  charge. 

2.2.2  Absorption  Results 

Before  discussing  the  emission  mechanism  and  results, 
it  is  useful  to  briefly  consider  how  below  bandgap  absorp¬ 
tion  might  occur  with  this  lineup.  Figure  20  shows  the 
p-In.  Ga  P/n-GaAs  hetero junction  in  equilibrium.  Notice 


that  this  configuration  yields  an  equilibrium  diagram 
very  different  from  that  of  p-ZnSnP2/n-GaAs .  In  the 
present  case  there  are  no  potential  wells  formed  adjacent 
to  the  interface.  Even  so,  a  possible  mechanism  for  below 
bandgap  absorption  still  can  be  seen.  An  electron  in  the 
valence  band  of  In.  Ga  P  may  absorb  a  photon  and  be 
excited  to  a  virtual  state  and  then  tunnel  into  the  con¬ 
duction  band  of  GaAs.  Another  possibility  would  be  an 
electron  tunneling  out  of  the  valence  band  of  In.  Ga  P  to 
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a  virtual  state  followed  by  photoexcitation  to  the  GaAs 
conduction  band.  The  latter  process  is  equivalent  to  hole 
tunneling  and  so  will  be  less  probable  due  to  the  larger 
hole  effective  mass.  Room  temperature  absorption  measure¬ 
ments  were  performed  on  these  heterostructures  in  the  same 

manner  as  for  the  ZnSnP-/GaAs  devices.  For  the  In.  Ga  P 
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on  GaAs  samples,  a  below  bandgap  peak  was  observed  in 
addition  to  the  peaks  associated  with  the  band-to-band 
absorption  in  the  two  materials.  The  below  bandgap  peak 
had  its  half-power  point  at  1.31  eV.  This  energy  cor¬ 
responded  exactly  with  the  predicted  value  of  Ej ,  however, 
this  exact  agreement  should  not  be  taken  too  seriously. 

The  below  bandgap  peak  increased  significantly  with 
reverse  bias.  This  is  expected  for  a  tunneling  assisted 
process  since  the  field  in  the  interfacial  region  increases 
with  increasing  reverse  bias. 


2.2.3  Emission  Results 

To  understand  how  this  hetero junction  might  act  as  an 
efficient  source  of  below  bandgap  emission,  look  at  the 
forward  bias  diagram  of  Figure  21.  Since  the  equilibrium 
lineup  contained  no  potential  wells  for  confinement,  these 
must  be  created  upon  the  application  of  sufficient  forward 
bias.  The  potential  wells  are  formed  by  the  accumulation 
of  majority  carriers.  With  this  n (down) /p (up)  lineup,  the 
below  gap  emission  mechanism  is  essentially  the  same  as  for 
the  p (down) /ri (up)  structure.  The  difference  is  that  accu¬ 
mulated  majority  carriers  instead  of  injected  minority 
carriers  are  the  species  involved  in  the  radiative  recom¬ 
bination.  Therefore,  the  emission  should  be  stronger 
since  there  is  a  more  plentiful  supply  of  recombining 
species. 

Room  temperature  emission  experiments  were  performed 
on  these  (p  on  n)  p-In.  Ga  P/n-GaAs  heterostructures.  The 
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samples  were  biased  at  constant  current  and  examined  with 

an  IR  viewer.  All  of  these  devices  yielded  little  or  no 

emission.  The  current  typically  reached  300  mA  or 
2 

150  A/ cm  before  any  emission  was  observed.  When  the 
#  o 

current  density  reached  about  250  A/cm  ,  the  I-V  charac¬ 
teristics  began  to  degrade,  probably  indicating  heating  of 
the  device.  In  all  cases,  the  emission  from  these  devices 
was  very  much  weaker  than  even  the  worst  of  the 


Figure  21  p-Ini-xGaxP/n-GaAs  forward  bias 

diagram  and  below  bandgap  emission 
mechanism. 


ZnSnP2/GaAs  devices.  For  comparison,  the  forward  bias 
emission  characteristics  of  (n  on  p)  n-In,  Ga  P/p-GaAs 
devices  were  also  studied.  Emission  was  observed  with 
the  IR  viewer  at  lower  current  bias  levels  than  for  the 
p  on  n  samples.  Although  the  emission  from  the  n  on  p 
samples  was  still  smaller  than  that  from  the  ZnSnP2/GaAs 
devices,  it  was  as  good  as  or  better  than  that  from  the 
p  on  n  devices. 

Our  results  show  that  for  the  In,  Ga  P/GaAs  mate- 
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rial's  system,  the  n (down) /p (up)  structure  did  not  prove 
to  be  the  strong  below  bandgap  emitter  as  predicted.  The 
emission  was  not  even  strong  enough  to  obtain  actual  spec¬ 
tra  through  the  monochromator.  The  emission  could  just 
barely  be  seen  with  the  IR  viewer.  Nothing  could  be  ob¬ 
served  in  this  manner  when  a  GaAs  filter  was  placed  in 
front  of  the  sample,  so  we  cannot  even  be  sure  that  below 
bandgap  emission  was  present.  We  believe  a  possible  reason 
that  this  configuration  did  not  prove  to  be  a  very  good 
emitter  is  linked  to  the  size  of  the  band  offsets.  To 
create  the  potential  wells  of  accumulated  majority  car¬ 
riers  on  either  side  of  the  interface,  a  large  forward 
bias  must  be  applied.  This  forward  bias  voltage  corre¬ 
sponds  to  an  energy  greater  than  either  of  the  bandgaps. 
This,  of  course,  is  not  possible  in  the  case  of  a  homo¬ 
junction,  because  this  would  entail  reaching  and  passing 


the  flat  band  condition.  At  flat  band,  for  the  homo¬ 
junction,  there  would  be  no  barrier  to  the  majority 
carrier  diffusion  and  so  the  junction  behaves  like  a 
pure  resistor.  The  homo junction  can  only  be  biased  up  to 
flat  band  because  beyond  this  point  the  heating  due  to 
the  high  current  would  destroy  the  device.  This  of  course 
is  the  reason  that  the  built-in  potential  is  the  maximum 
voltage  that  can  be  applied  safely  to  a  homo junction.  If 
there  are  barriers  to  majority  carrier  diffusion  in  the 
form  of  band  offsets,  this  flat  band  condition  in  theory 
could  be  reached  and  passed.  Recall  that  our  measured 
band  offsets  for  the  p-In.  Ga  P/n-GaAs  devices  were 
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AEC  =  0.59  eV  and  AEV  =  0.08  eV.  Despite  the  fact  that 
these  values  may  not  be  very  accurate,  the  valence  band 
offset  appears  to  be  small.  Below  some  critical  value  of 
either  band  offset,  which  is  related  to  doping  and  other 
material  parameters,  prohibitively  large  currents  will 
flow  in  forward  bias.  The  hetero junction  in  this  case  will 
not  be  much  different  from  a  homo junction.  Therefore, 
one  will  not  be  able  to  go  sufficiently  beyond  the  flat 
band  condition  to  create  the  charge  accumulation  wells. 

This  is  probably  the  reason  our  p  on  n  devices  proved 
to  be  such  poor  emitters.  Also  the  n  on  p  emitters  may 
have  been  generally  inferior  to  those  made  of  ZnSnP2/GaAs 
since  the  carrier  concentrations  in  the  latter  case  were 
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much  higher.  The  electroluminescence  observed  by  Caine 
et  al.  [3]  from  the  p-Al  In,  As/n-InP  devices  at  1.4K 
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indicated  that  these  were  better  emitters  than  our 
n (down) /p (up)  devices.  The  band  offsets  for  these 
p-Al  In..  As/n-InP  hetero junctions  were  determined  to 
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be  AE  =  0.52  eV  and  AE  =  0.40  eV  .  Despite  the  inaccu- 
c  v 

racy  of  band  offset  measurements,  these  values  may  indi¬ 
cate  large  enough  barriers  to  current  flow  to  allow  these 
devices  to  be  biased  beyond  the  flat  band  condition. 
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2.3  C-V  INTERCEPT  METHOD 


Theoretical  models  for  the  capacitance  of  an  abrupt 
p-n  hetero junction  can  be  made  very  complex.  As  men¬ 
tioned  earlier,  the  expression  used  to  interpret  the  C-V 
data  for  the  ZnSnP2/GaAs  heterostructures  neglected  the 
effects  of  interfacial  charge  and  dipoles.  That  expres¬ 
sion,  equation  (2.3),  is  valid  for  an  abrupt  hetero junction 
in  the  depletion  approximation.  To  modify  this 
model  of  the  capacitance,  let  us  consider  an  abrupt  hetero¬ 
junction  with  a  sheet  of  positive  charge  at  the  interface. 
For  interface  charge  density  qa5(0),  Donnelly  and 
Milnes  [18]  obtain. 


r- 2  _  2 (eaNa  +  gnNd) (VQ  -  V)  -  qg2  (4.1) 

qenepAzNdNa  * 
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This  is  the  same  expression  for  C  as  quoted  earlier  in 
equation  (2.3),  except  for  an  additional  negative  term  in 
the  numerator.  Equation  (4.1)  reveals  that  our  previous 
C-V  intercept  determination  of  VQ  neglecting  interface 

charge  will  be  in  error  if  significant  interface  charge 
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is  present.  From  (4.1)  we  see  that  C  =0  when. 
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Here  V  is  the  true  built-in  potential  whereas  V  is  the 


measured  value.  This  expression  reveals  that  regardless 
of  the  sign  of  any  interface  charge,  the  measured 
built-in  potential  will  be  smaller  than  the  real  value. 

For  this  reason,  we  used  the  largest  value  of  V0m,  pre¬ 
sumably  the  most  accurate  value  for  constructing  the 
p-2nSnP2/n-GaAs  band  lineup.  As  can  be  seen  from  (4.1), 

this  voltage  independent  interfacial  charge  will  not 
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affect  the  slope  of  the  C  -V  plot  and  so  linearity  of 

this  plot  will  not  guarantee  that  interface  charge  is  not 
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a  problem.  A  nonlinear  C  -V  plot  may,  however,  indicate 

several  other  properties  of  the  interface.  Two  inter facia 

effects  [18]  which  would  lead  to  a  nonlinear 
-2 

C  -V  plot  are  voltage  dependent  interfacial  charge 
and  dipoles. 

A  voltage  independent  interfacial  dipole  would  not 

_2 

affect  the  slope  of  the  C  -V  plot  but  would  change  the 
intercept  and  so  VQm.  VQm  could  be  decreased  or  increased 
from  the  dipole  free  value  depending  on  the  sign  of  the 
dipole.  A  positive  interfacial  dipole  is  defined  as  that 
which  serves  to  increase  the  measured  value  V0m.  This 
constant  dipole  may  vary  from  sample  to  sample.  Unfor¬ 
tunately  this  dipole  changes  not  only  V0m,  but  also  VQ 
and  so  the  band  offsets.  For  simplicity  we  will  not  con¬ 
cern  ourselves  further  here  with  constant  interfacial 
dipoles . 
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Voltage  dependent  interfacial  dipoles  and  charge  will 
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contribute  to  a  variable  slope  of  the  C  -V  plot.  Depend¬ 
ing  on  the  rate  of  change  of  the  dipole  and  charge  with 

_2 

voltage,  the  C  -V  plot  will  curve  upward  or  downward 

~ith  increasing  reverse  bias.  The  exact  details  of  the 
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behavior  of  the  C  -V  plot  in  the  presence  of  voltage 
dependent  interfacial  dipoles  and  charge  will  not  be 
dealt  with  here,  but  can  be  found  in  a  previous  work  [18] 

We  have  been  discussing  the  effects  of  interfacial 
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charge  and  dipoles  on  the  C  -V  plot  of  an  abrupt  hetero¬ 
junction.  Now  let  us  examine  the  characteristics  of 
these  plots  for  the  ZnSnP,/GaAs  and  In1  Ga  P/GaAs  devices 

used  in  our  below  bandgap  optical  measurements. 

-2 

Firstly,  all  of  the  1  MHz  C  -V  plots  for  the 
ZnSnP2/GaAs  samples  studied  were  straight  lines  with 
nearly  the  same  slope.  The  built-in  potential  estimates 
from  the  intercepts  of  these  lines  ranged  from  1.33  to 
1.75  V.  As  already  mentioned,  all  V0jn  values  except  for 
V0m=  1.33  V  confirm  a  staggered  lineup  for  this  hetero¬ 
junction.  Though  I-V  characteristics  are  not  generally 
very  precise  measures  of  interfacial  properties,  for 
certain  hetero junctions  the  I-V  characteristics  can 
still  indicate  a  lower  limit  on  V  .  The  I-V  curve,  for 
the  sample  yielding  VQm=  1.33  V,  revealsd  a  Vq  at  least 


greater  than  1.6  V.  Guided  by  this  result,  let  us  suppose 

for  simplicity  that  the  variability  in  the  values  is 

due  to  interface  charge.  Assuming  that  VQ=  1.75  V,  we  can 

calculate  the  interface  charge  density  necessary  to  shift 

V0m  to  1.33  V.  If  one  inserts  the  appropriate  values  for 

the  materials  parameters  quoted  earlier  into  (4.2),  we 
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obtain  the  interface  charge  density,  a  =  1.5  x  10  cm 

Though  several  mechanisms  could  create  an  interface 
charge  of  this  magnitude,  let  us  consider  the  charge  that 
would  result  for  slight  misorientations  from  a  charge 
neutral  orientation.  Earlier,  an  outline  of  Harrison ' s  [ 7  ] 
arguments  concerning  charge  neutral  hetero junction  inter¬ 
faces  was  given.  The  {211}  orientation  used  for  the  GaAs 
substrate  in  the  ZnSnP2/GaAs  hetero junctions  was  shown  to 
be  theoretically  charge  neutral.  Using  Harrison's  simpli¬ 
fied  model,  even  slight  misorientations  from  the  charge 
neutral  orientation  would  result  in  considerable  interface 
charge.  It  can  be  shown  that  the  charge  density  that 
results  for  small  misorientations  from  {211}  is, 

qa  =  (q/T/a2)sinO  .  (4.3) 

Here  a  is  the  lattice  constant  and  0  is  the  angle  of  mis- 
orientation  [20],  The  practical  limit  of  current  routine 
wafer  orientation  techniques  is  roughly  ±0.5°.  For  0.5° 
misorientation ,  using  a  =  5.6532  A  for  GaAs,  we  obtain 
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a  =  4.7  x  10  cm  .  Since  this  charge  is  not  enough  to 

account  for  a  measurement  error  in  V  of  0.42  V,  there 

o 

may  be  other  interfacial  charging  mechanisms  involved  in 
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these  devices.  Next,  let  us  consider  the  C  -V  plots  for 
the  In.  Ga  P/GaAs  devices.  Both  lattice-matched  and 

X 

mismatched  samples  were  investigated.  These  devices  were 
all  {lll}p-In.  Ga  P/n-GaAs.  At  1  MHz,  the  C~2-V  plots 
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for  the  lattice-matched  devices  were  straight  lines  with 

some  variability  in  the  intercept.  The  mismatched  sam- 
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pies  had  C  -V  plots  at  1  MHz  which  curved  downwards  for 
increasing  reverse  bias.  These  results  indicate  that 
the  mismatch  in  these  devices  may  create  voltage-dependent 
interfacial  dipoles  and  charges.  The  lattice  mismatch  in 
these  samples  was  approximately  0.5%  which  is  only  slightly 
greater  than  the  0.44%  mismatch  associated  with  the 
ZnSnP2/GaAs  devices.  The  {111}  interface  of  these 
In.  Ga  P/GaAs  devices  is  not,  however,  theoretically  charge 

a  X 

neutral  because  of  electronegativity  differences.  Therefore, 
the  mismatch  and  the  orientation  together  may  serve  to  create 
more  significant  interface  charge  and  dipole  effects  than 
found  in  the  ZnSnP2/GaAs  devices. 

_2 

Another  interesting  effect  was  seen  in  C  -V  plots 
for  measurements  taken  at  1  kHz.  Initially  C-V  measure¬ 
ments  were  performed  at  1  kHz  on  the  ZnSnP2/GaAs  devices 


J 
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_2 

using  the  HP  LCR  meter.  These  C-V  data  yielded  a  C  -V 

plot  with  a  decreasing  slope  at  increasing  reverse  bias, 

shown  in  Figure  22.  The  extrapolation  of  the  reasonably 

linear  portion  of  this  plot  at  low  bias  yielded  a  higher 

value  of  than  for  the  same  sample  at  1  MHz.  The 

problem  may  lie  within  the  measurement  technique  and  not 

the  sample.  The  HP  meter  uses  a  50mV  rms  test  voltage. 

This  may  be  too  large  for  small  signal  approximations  to 

be  valid.  When  these  same  devices  were  measured  at  1  kHz 

with  a  15mV  rms  test  voltage  using  the  lock-in  technique, 

-2 

linear  C  -V  plots  were  obtained.  A  typical  result  is 

shown  as  the  solid  line  in  Figure  22.  The  voltage  inter- 

_  9 

cept  from  this  C  -V  plot  at  1  kHz  yielded  about  the  same 
value  of  as  for  this  sample  at  1  MHz. 

Apparent  from  these  results  is  the  great  care  needed 
to  perform  meaningful  capacitance  measurements  on  hetero¬ 
junctions  and  also  the  difficulty  in  interpreting  these 
results.  Various  types  of  capacitance  measurements  will 
undoubtedly  continue  to  be  developed  as  tools  for  analyzing 
hetero junction  properties. 
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2 . 4  CONCLUSIONS 

This  work  has  sought  to  investigate  the  predicted 
phenomena  of  below  bandgap  emission  and  absorption  in 
staggered  lineup  hetero junctions .  These  phenomena  were 
predicted  based  upon  a  model  involving  mechanisms  pecu¬ 
liar  to  the  staggered  lineup.  We  observed  below  bandgap 
absorption  in  two  materials  systems,  ZnSnP2/GaAs  and 
In^^Ga^/GaAs  and  below  bandgap  emission  in  at  least 
the  first  of  these  systems.  Another  research  group  [4] 
also  observed  below  bandgap  electroluminescence  at  low 
tempera ture  in  Al.  In  As/InP.  Instead  of  reviewing 

X  ™  A  X 

these  results  in  detail,  we  will  discuss  some  of 
the  difficulties  in  our  interpretations  of  these  data. 

We  have  already  discussed  some  of  the  pitfalls 
in  the  experimental  determination  of  band  offsets  in 
heterostructures.  One  difficulty  not  previously 
mentioned  is  that  of  actually  growing  abrupt 
hetero junctions.  Our  results  were  all  taken  from 
samples  grown  by  liquid  phase  epitaxy.  There  is 
skepticism  about  the  ability  of  this  growth  technique 
to  achieve  sufficiently  abrupt  interfaces  for 
device  applications  such  as  ours.  There  are 
other  growth  techniques,  such  as  molecular  beam 
epitaxy  (MBE)  and  organo-metallic  vapor 


phase  epitaxy  (OMVPE) ,  which  have  achieved  more  abrupt 
hetero-interfaces.  We  do  not  know  how  abrupt  our  hetero¬ 
junctions  are  and  do  not  know  even  how  critical  abruptness 
might  be  to  the  mechanisms  we  have  described.  That  is,  the 
devices  might  still  work  if  the  interfaces  were  significantly 
graded,  but  the  analysis  would  be  much  more  complicated. 

Another  difficulty  with  our  interpretation  is  the 
possibility  that  tunneling  assisted  transitions  involving 
impurity  or  defect  levels  at  the  interface  could  account 
for  our  results.  This  alternate  explanation  addresses  the 
particular  problem  of  the  shape  of  the  below  bandgap  photo¬ 
response.  In  both  the  emission  and  absorption  spectra,  the 
below  bandgap  transition  manifests  itself  in  the  form  of  a 
peak.  One  can  easily  account  for  the  peak  in  the  emission 
spectrum.  Emission  involves  transitions  between  filled 
states  at  higher  electron  energies  and  unfilled  states  at 
lower  energies.  When  considering  the  initial  filled  states 
in  the  conduction  band  well,  one  sees  that  these  states  are 
less  likely  to  be  occupied  at  higher  energies.  Also,  the 
states  in  the  valence  band  well  are  more  likely  to  be 
filled  the  further  down  in  energy  one  looks  into  the 
valence  band.  Both  these  effects  give  qualitative  plausi¬ 
bility  for  a  peak  rather  than  a  shoulder  appearing  in  the 
below  bandgap  emission. 


The  below  bandgap  absorption  peak,  however,  is  not 
so  easy  to  explain.  In  this  case,  one  starts  with  a  filled 
state  in  the  valence  band  well  and  an  empty  state  in  the 
conduction  band  well.  All  of  the  obvious  factors  govern¬ 
ing  this  type  of  transition;  occupation  numbers,  density 
of  states,  and  tunneling  probabilities  serve  to  increase 
the  response  at  higher  photon  energies.  There  is,  there¬ 
fore,  no  obvious  mechanism  that  would  cut  down  the  response 
between  the  minimum  energy  Ej  and  the  smaller  of  the  two 
bandgap  energies.  We  would,  therefore,  expect  to  see  a 
shoulder  in  the  photoresponse  instead  of  a  peak. 

A  possible  alternate  explanation  of  our  results 
might  be  tunneling  assisted  optical  transitions  involving 
impurity  or  defect  levels  at  or  near  the  interface.  The  density  of 
states  from  these  may  peak  in  some  region  of  the  gap  and 
could  easily  explain  both  the  below  bandgap  absorption 
and  emission  peaks.  We  looked  for  evidence  in  our  data 
contrary  to  this  interpretation.  We  observed  that  the 
below  bandgap  emission  and  absorption  peaks  were  coinci¬ 
dent  in  energy.  In  principle  one  expects  a  shift  between 
emission  and  absorption  spectra  for  transitions  involving 
impurities  because  of  the  Franck-Condon  principle.  As 
discussed  earlier,  this  principle  states  that  in  optical 
transitions,  electrons  move  so  quickly  that  the  slower 


heavier  lattice  atoms  can  be  considered  static  during  the 
transition.  The  lattice  relaxation  tends  to  result  in 
smaller  energies  in  emission  for  such  transitions  than  in 
absorption.  Franck-Condon  shifts  between  emission  and 
absorption  spectra  for  different  impurity  levels  in  various 
semiconductors  have  been  measured.  Unfortunately,  no  re¬ 
liable  data  are  available  that  could  be  directly  used  as  a 
guide  to  estimate  the  shift  in  our  case.  It  is  possible 
that  this  shift  may  be  too  small  to  be  observed  considering 
the  broadness  of  our  spectral  peaks.  Another  possible 
cause  of  the  below  bandgap  peak  in  the  absorption  may  be 
exciton  absorption  at  the  interface.  Exciton  levels  can 
give  rise  to  a  peak  in  the  density  of  states  in  the  gap. 
Room  temperature  absorption  measurements  of  various  semi¬ 
conductors  have  revealed  below  bandgap  peaks  due  to  exci- 
tonic  transitions. 

These  are  some  of  the  problem  areas  in  our  investi¬ 
gation  of  below  bandgap  emission  and  absorption  in 
staggered  lineup  heteroj unctions .  Although  there  is 
evidence  that  this  type  of  heteroj unction  may  prove 
useful  through  the  mechanisms  we  have  described,  con¬ 
siderably  more  research  into  this  area  is  needed  before 
any  solid  conclusions  could  be  reached. 
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3.  SPATIALLY  VARYING  ENERGY  BANDS 


With  the  development  of  special  growth  techniques  such  as  molecular  beam 
epitaxy  and  metal  organic  chemical  vapour  deposition,  tailor-made  submicron 
semiconductor  heterostructures  can  now  be  designed.  An  entire  new  variety  of 
structures  including  quantum  wells,  superlattices  and  modulation-doped  struc¬ 
tures  can  be  conceived,  and  some  have  already  led  to  successful  device  appli¬ 
cations. 

To  quote  a  few,  quantum  wells  are  being  investigated  in  semiconductor 
lasers  [1]  and  resonant  tunneling  [2]  for  submillimeter  radiation.  Super- 
lattices'  potential  applications  Include  non-linear  optics  in  relation  with 
optical  switching  (bistability)  [3],  the  Bloch  oscillator  [4]  and  the  Zener 
oscillator  [5]  for  submillimeter  radiation.  An  Important  application  of  modu¬ 
lation  doping  [6]  is  the  high  electron  mobility  transistor,  a  low-noise  micro- 
wave  amplifier  [7]  and  a  candidate  for  high  speed  logic  [8]. 

These  semiconductor  heterostructures  are  made  of  successive  semiconductor 
crystal  layers  grown  on  top  of  each  other.  Provided  the  lattice  parameters 
are  closely  matched,  the  lattices  of  the  semiconductors  essentially  cohere 
with  a  minimal  perturbation.  Due  to  the  special  growth  techniques  used,  the 
spatial  variation  of  the  semiconductor  materials  and/or  of  the  doping  can  be 
controlled  so  as  to  occur  in  a  few  lattice  parameters. 

The  study  of  the  physical  properties  of  these  semiconductors'  structures 
is  the  object  of  extensive  theoretical  and  experimental  effort  in  relation 
with  the  important  potential  device  applications,  some  of  which  were  quoted 
above,  most  of  which  remain  however  to  be  devised. 

A  fundamental  property  on  which  all  these  heterostructures  rely  is  the 
spatial  variation  of  the  bandgap  and  the  conduction  and  valence  bandstruc- 
tures.  In  this  work  we  shall  be  concerned  with  the  development  of  a  consis¬ 
tent  and  simplified  picture  of  spatially  varying  bandstructures. 


Bandstructures  derive  from  the  periodic  nature  of  the  crystal  potential 


and  are  therefore  mathematically  well-defined  only  for  an  Infinite  crystal. 
However,  In  practice,  the  concept  of  bands  holds  on  a  microscopic  scale. 

This  Is  exemplified  by  the  use  of  band  diagrams  In  classical  device  theory. 

The  total  classical  Hamiltonian  H  can  Indeed  be  written  as 

H(i?.r)  -  E(ii)  -  eV(?)  (1) 

where  E(£)  Is  the  conduction  or  valence  bandstructure,  e  the  electron  charge 
and  V(r)  the  electrostatic  potential.  The  symbol  r  denotes  the  spatial 
location  and  £  the  Bloch  wavevector  or  quasi -momentum.  In  a  band  diagram 
usually  only  the  bottom  of  the  conduction  band  and  the  top  of  the  valence  band 
are  represented,  which  Is  equivalent  to  specifying  £. 

Equation  (1)  together  with  the  acceleration  theorem  [9]  constitute  the 
semi -classical  picture  on  the  basis  of  which  one  can  describe  the  ballistic 
motion  of  an  electron  or  derive  (using  the  Boltzman  Transport  Equation)  the 
standard  device  equations  [10].  In  the  heterostructures  discussed  above,  the 
spatial  dimensions  are  made  sufficiently  small  for  the  quantization  of  the 
electron  states  to  be  significant.  One  needs  then  to  revert  to  the  quantum 
mechanical  form  of  equation  (1)  where  the  quasi -momentum  Is  replaced  by  the 
operator 

k  *  -1  ~  (2) 

3r 

«• 

and  to  solve  the  Schrodlnger  equation 


ttl  ■57  *(t,f)  ■  H  *(t,f) 


Equations  (1)  together  with  (2)  are  often  referred  to  as  the  Wannier 
theorem,  the  derivation  of  which  can  be  found  in  textbooks  on  semiconductor 
theory  [11],  The  Wannier  theorem  relies  on  two  basic  assumptions.  The  first 
is  the  one-band  approximation  which  holds  If  the  perturbation  potential  V(r) 
is  not  strong  enough  for  band-mixing  to  occur  or  Interband  transition  such  as 
tunneling  to  take  place.  Secondly,  the  perturbation  potential  should  vary 
slowly  or  smoothly  with  position  so  as  to  constrain  the  potential  overlap 
matrix  to  a  diagonal  representation. 

These  assumptions  are  adequate  for  the  classical  representation  in  which 
the  electrons  are  essentially  free  and  can  be  represented  by  a  wavepacket  j(r). 
For  the  submicron  heterostructures  of  Interest  here,  these  assumptions  do  not 
hold,  as  the  potential  variations  In  modulation  doped  or  in  spatially  varying 
bandgap  structures  are  both  large  and  can  occur  in  a  few  lattice  parameters. 

In  this  paper  we  shall  propose  a  concept  of  generalized  bandstructures 
for  which  the  Wannier  theorem  can  be  generalized  and  applied  to  submicron 
heterostructures. 

First,  we  shall  rederive  the  Wannier  theorem  in  Section  3.1  specializing  it 
to  one-dimensional  structures.  Since  the  submicron  heterostructures  grown  by 
MBE  or  MOCVD  are  layered  structures,  a  one-dimensional  picture  can  be  imple¬ 
mented  provided  some  assumptions  are  made.  Such  one-dimensional  models  are  of 
heuristic  interest  as  their  simplicity  promotes  better  insight,  and  they  are 
sometimes  used  in  semi -quantitative  studies. 

In  Section  3.2,  we  discuss  the  extension  of  the  Wannier  theorem  using  the 
Generalized  Wannier  Functions  and  study  in  Section  3.3  two  examples  of  hetero- 
structures.  These  examples  constitute  both  a  test  and  a  demonstration  of 
applications  of  the  generalized  Wannier  picture. 
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Final  ly,  we  present  in  Section  3.4  the  application  of  the  generalized 
Wannler  picture  to  "true”  heterostructures  and  discuss  the  generalization  of 
the  concept  of  bands  to  spatially  varying  bandstructures: 

3.1  ONE-DIMENSIONAL  WANNIER  MODEL 

In  this  section,  we  consider  the  problem  of  a  crystal  electron  In  an 
external  one-dimensional  potential  U.  The  total  Hamiltonian  of  the  electron 
Is  given  by: 

H  -  H0  +  U(d«r) 

where  H0  is  the  unperturbed  crystal  Hamiltonian  and  U(x)  Is  the 
one-dimensional  energy  potential  varying  along  an  axis  (the  device  axis) 

A 

represented  by  the  unit  vector  d. 

We  assume  the  device  axis  to  be  parallel  to  the  lattice  vector.  We  only 

consider  face-centered  cubic  semiconductor  lattices  for  which  a  lattice  vector 
♦  ♦ 

R(n)  can  be  written  using  the  orthonormal  basis  of  the  Bravals  lattice  as 
*  ♦  ♦  ♦ 

R(n)  *  An  with  n  a  set  of  Integers  (nlt  n2,  n3),  a*  the  lattic  parameter 
of  the  cubic  lattice  and  A  the  matrix 


♦ 

We  shall  denote  as  a  the  modulus  of  the  smallest  lattice  vector  a  parallel  to 

A  A 

d;  a  Is  therefore  the  effective  lattice  parameter  along  the  direction  d  and 

*  ♦  ♦ 

any  lattice  vector  parallel  to  d  is  written  R  »  na  with  n  an  integer. 

♦  ♦  ♦  ♦ 
A  reciprocal  lattice  vector  X  is  written  in  the  same  basis  as  K(l)  *  B1 

with  1  a  set  of  integers  (1^  12,  13)  and  B  the  matrix 


"J  '  v 


:r  ■  v*'  v  >■ ’.’ 


L"  I  --  -■  ’VT  LP  "-F  ’ 
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So  that  AB  *  2*.  It  Is  easy  to  prove  that  for  cubic  lattices,  there  always 

exists  a  reciprocal  lattice  vector  parallel  to  a  given  direct  lattice  vector. 

♦  ♦ 

This  Is  equivalent  to  finding  the  constant  c  such  tht  K  *  c  R.  A  possible 
solution  Is  c  ■  8*/a*2#  which  leads  to 


r2 


: 


r  r 

2  1 

1  2. 


n. 


♦ 

We  shall  denote  q  the  modulus  of  the  smallest  reciprocal  lattice  vector  q 
*  * 

parallel  to  d  so  that  any  reciprocal  vector  parallel  to  d  can  be  written 
♦  + 

K  *  lq  with  1  an  Integer.  Since  the  lattice  vector  and  the  reciprocal  lattice 
vector  satisfy  the  relation 


♦  ♦  ♦  ♦  ♦  ♦ 
R(n)  •  K(1 )  *  2*  n  •  1 


We  have  that  q  and  a  satisfy 
♦  ♦ 

a  .  q  •  aq  »  2irp 


«►  ♦ 

with  p  an  Integer.  For  the  <100>  direction  (R  ■  a( 1,0,0) )  we  have  n  *  (-1,1,1) 
and  the  smallest  reciprocal  vector  is  given  by  1  *  (0,1,1),  so  that  a»a*  and 
P-2. 

♦  ♦ 

In  the  extended  zone  scheme  the  Bloch  function  solution  *(k,r)  of 

(3) 


H0  ♦0c)  •  E(k)  *(k) 


Is  a  periodic  function  of  k  along  the  device  direction  d.  We  can  introduce  a 

♦  a 

one-dimensional  quasi  momentum  k  along  the  device  direction  d  defined  by 
k(k)  «  kx  +  k  d 

♦  « 
with  kxthe  transverse  momentum  perpendicular  to  the  direction  d.  It  follows 

A 

that  a  Bloch  function  along  the  direction  d  can  be  written  as 
*  +  + 

♦(kd  +  kitr)  »  $(k)  which  is  a  periodic  function  of  k  with  period  q: 

♦(k  +  nq)  *  $(k)  for  n  an  integer. 

A 

We  can  now  define,  along  the  direction  d,  a  one-dimensional  Wannier  Func¬ 
tion  w(m,r,icx)  as  the  spatial  Fourier  coefficient  of  the  Bloch  function  $(k) 

r-  q/Z  -  *  ♦ 

w(m,r,kx)  *  l//q  /  ♦(kd  +  kx,r)  exp  (-1km2ir/q)  dk  (4) 

-q/2 

and  reciprocally  we  have 

♦  (kd  +  icx.r)  ■  l//q  £  w(m,r,iTx)  exp  (ikm2ir/q) 

m 

From  the  orthogonality  property  of  the  Bloch  functions  and  from  equation  (3) 
one  easily  verifies  that  these  Wannier  functions  (WF's)  form  an  orthogonal  set 

m  *  *  ■¥  ♦ 

<w(n)|w(m)>  »  /  w(n,r,kx)*  w(m,r,kx)  dx'  *  5nm 

with  x1  =  d*  r  and  satisfy  the  useful  relation  H  w(m,r  )  =  rE(n-m)w(n,r,k,) 
q/2  0  1  n  1 

with  E(n)  =  1/q  f  +  kd)  exp  (-i kn2it/q)  dk  . 
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The  E(n)  are  recognized  as  the  Fourier  coefficients  of  the  bandstructure  along 
♦ 

the  line  k(k)  parallel  to  the  device  axis. 

For  a  sufficiently  smooth  and  weak  one-dimensional  potential  U(d»r)  no 
band  mixing  Is  expected  and  the  transverse  momentum  Is  conserved.  One  can 
then  expand  the  electron  state  In  Wannler  Functions  of  the  band  considered 
with  an  envelope  f(n,t) 

J  ■  I  f(n,t)  w(n,r,ici)  (5) 

n 

The  weight  function  |f(n,t)|2  Is  the  probability  of  finding  the  electron  at 
the  lattice  site  n  at  time  t.  Replacing  t  In  the  Schrodlnger  equation 

Ht  *  **  It  T  (6) 

by  the  expansion  of  equation  (5).  multiplying  equation  (6)  by  w(n)*  and 
Integrating  over  x1,  one  obtains  the  Wannler  recurrence  equation 

l  [E(n-m)  ♦  U(n,m)]  f(m,t)  »  tfi  f(n,t)  (7) 

m 

where  U(n,m)  Is  the  matrix  element  given  by 

m 

U(n,m)  *  <w(n) |U|w(m)>  «  /  w{n,r,iTji)*  U(£»r)  w(mtrtiTi)  dx'  (8) 


,*-\v \*-\\ 
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The  one-band  Wannler  recurrence  equation  derived  holds  for  sufficiently  small 
perturbation  potentials  U.  For  equation  (7)  to  be  solved,  both  the  band- 
structure  Fourier  coefficients  and  the  matrix  elements  U(n,m)  are  required. 

As  expressed  by  equation  (8)  the  matrix  elements  are  evaluated  from  the 
Wannler  functions.  Techniques  for  the  evaluation  of  the  Wannier  functions 
have  been  developed  by  Kohn  [12],  The  evaluation  of  the  WF  for  silicon  was 
recently  reported  by  Evan  Kane  and  Anne  Kane  [13],  however,  the  WF  derived  is 
a  linear  superposition  of  the  WF's  of  the  four  valence  bands.  The  evaluation 
of  realistic  WF's  Is  a  complicated  matter  which  we  would  like  to  avoid  in  the 
simple  picture  we  Intend  to  develop. 

Some  direct  approximations  can  be  derived  from  the  properties  of  the 
WF's.  In  the  case  of  simple  bands  It  is  possible  to  select  the  phase  of  the 
Bloch  functions  so  that  the  WF's  are  exponentially  localized  [14],  From  both 
the  orthogonality  and  the  tight  localization  of  the  WF's,  it  follows  that  the 
matrix  element  of  a  smooth  potential  is  accurately  represented  by  the  sampling 
of  the  potential  at  the  lattice  sites 

U(n,m)  *  U(x'  *  na/p)  5n  m  (9) 

We  shal  1  discuss  in  Section  3. 3  the  extent  to  which  these  approximations 
can  be  used  for  sharply  varying  potentials. 

Finally,  an  Important  property  of  the  one-dimensional  WF  we  did  not 
address  so  far  Is  their  Invariance  under  translation  from  one  lattice  site  to 
another.  Following  the  approach  of  Wannier  [15],  we  first  write 

?(£,?-$)  -  l//q  l  exp[ikma/p]  wCm.r-tf.icx)  (10) 

m 

which  using  the  translation  property  of  the  Bloch  function  can  be  written 


T(i?,p-$)  *  l//q  l  exp[-1it«ff]  exp[1kma/p]  w(mtr,itjj 
m 

Fop  ^  *  n  a,  we  have  Ic  •  -  (kd  +  (T^)  •  J  *  k  d  •  I  *  knat  and  changing  the 

Index  of  summation  we  obtain 

T(it,r-ff)  *  l//q  l  exp[ikma/p]  w(m+np,r,iti)  (11) 

m 

Identifying  the  coefficients  in  the  series  (10)  and  (11)  we  finally  get 
w(m,r-na,iti)  «  w(nH-nptr,icj.). 

Considering  the  <100>  direction  for  which  p*2,  it  follows  that  there  exist  two 
sets  of  WF's  each  of  which  can  be  generated  from  two  generic  WF's  by  a  lattice 
translation 

w(2n,f)  *  w(0,r-na) 
w( l+2n,r)  ■  w(l,?-nl) 

The  number  of  WF's  for  a  given  direction  Is  therefore  given  by  p  for  the  face- 
centered  cubic  lattice  considered.  For  our  purpose  it  is  convenient  to 
picture  these  WF's  as  located  at  p  sites  along  the  lattice  vector  a  associated 

A 

with  the  direction  d.  We  can  then  write  these  WF's  as 
w(n,r,ici)  *  w(x'-na/p,ri,icjL)  *  w(x'-na/p)  »  w(n) 


where  a  Is  the  effective  lattice  parameter  for  the  direction  d.  The  one- 
dimensional  WF  we  have  introduced  is  in  fact  a  hybrid  entity,  since  it 


corresponds  to  the  Bloch  State  k^  for  the  transverse  coordinates  rj_.  The 
properties  of  the  WF’s  of  one-dimensional  lattices  cannot  be  assumed  to  apply 
directly  to  the  one-dimensional  WF  we  have  defined.  Therefore,  problems 
concerning  the  reality  of  the  WF  and  the  simultaneous  convergence  of  the 
WF's  remain  to  be  addressed.  These  studies  are  not  trivial  matters  for 
three  dimensional  crystals  [14],  In  particular  for  lattices  without  a  center 
of  Inversion. 

In  our  simplified  treatment  we  shall  assume  the  picture  developed  to  be 
truly  one-dimensional.  The  reality  and  exponential  localization  of  the  WF's 
in  one-dimensional  structures  are  discussed  by  Kohn  [16].  More  recent  work  by 
Zak  [17]  raises  the  one-dimensional  Wannier  picture  to  the  status  of  a  well- 
defined  representation  In  quantum  mechanics  and  introduces  the  canonical 
Wannier  function  with  minimal  position  uncertainty  located  in  the  zero  cell 
of  the  Bravais  lattice. 

Zak  also  demonstrates  that  the  position  operator  has  a  discrete  spectrum. 
This  implies  that  the  applied  potentials  are  transformed  by  the  one-band 
approximation  into  step-like  potentials.  This  is  due  to  the  inability  of  the 
one-band  approximation  to  deal  with  the  perturbation  of  the  inner  crystal 
potential.  However,  as  we  shall  see  in  the  next  section,  most  of  the  band 
mixing  can  be  formally  accounted  for,  with  the  use  of  the  generalized  Wannier 
function.  For  stronger  potentials  it  is  further  necessary  to  account  for  the 
polarization  of  the  semiconductor  [18], 
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3.2  GENERAL  WANNIER  MODEL 

In  this  section  we  shall  be  concerned  only  with  abrupt  modulation -doped 
heterostructures  or  models  of  abrupt  heterostructures  for  which  the  bandgap  Is 
spatially  varying  while  the  bandstructure  Is  assumed  to  be  conserved. 

The  one-dimensional  Wannler  picture  developed  In  the  previous  section  Is 
capable  of  handling  fast  spatial  variation  of  the  potential  since  the  non- 
diagonal  term  of  the  matrix  potential  V(n,m)  can  be  accounted  for.  However, 
the  formalism  developed  relies  on  the  one-band  approximation  which  strictly 
holds  only  for  small  potential  perturbations.  Its  use  for  submicron  hetero- 
structures  seems  unjustifiable.  Indeed,  the  band-edge  variation  In  these 
structures  Is  large  and  band  mixing  Is  expected. 

One  should  then  resort  to  a  multiband  expansion,  and  the  envelope 
equation  (7)  now  reads 

til  f(n,b)  -  l  {E(n-m,b)  f(m,b)  +  J  <w(n,b) |U(x) |w(m,b' )>  f(m,b')} 
m  b 1 

where  b  is  the  band  Index. 

It  Is,  however,  possible  to  avoid  a  multiband  expansion  If  one  Introduces 
the  concept  of  generalized  bands.  As  was  remarked  by  Zak  [17],  "a  band  Is  a 
surprisingly  stable  entity  under  external  perturbations."  This  might  seem  to 
be  a  paradoxical  statement  coming  from  Zak,  who  has  been  very  critical  of  the 
use  of  the  one-band  approximation  In  general  [19],  and  in  particular  for  the 
Stark  effect  [20],  In  fact,  the  existence  of  bands  under  external  perturba¬ 
tions  does  not  rely  on  the  one-band  approximations.  A  band  of  energy  Is 
essentially  a  quasi -continuum  of  energy  states.  Such  an  entity  seems  to  sur¬ 
vive  numerous  mistreatments  such  as  high  impurity  concentrations,  alloying 
with  random  distribution  of  two  atoms  on  a  given  sublattice  In  III-V  compounds 


and  even  abrupt  hetero junctions.  However,  surface-states  or  band  tailing  or 
even  alloy  and  roughness  scattering  effects  occur*  and  It  might  be  more  accu¬ 
rate  to  refer  to  these  bands  as  generalized  bands.  Such  bands  In  most  cases 
can  only  be  represented  mathematically  as  resulting  from  the  mixing  of  the 
bands  of  the  unperturbed  crystals.  This  Is  similar  to  the  expansion  of  the 
pseudo-wavefunctlon  In  terms  of  plane  waves  In  the  calculation  of  band- 
structures  using  the  pseudo-potential  method;  the  bands  evaluated  result  from 
the  mixing  of  the  bands  of  the  empty  lattice. 

A  one-band  formalism  can  then  be  developed  If  one  uses  the  exact  Wannler 
functions  instead  of  the  Wannler  functions  of  the  unperturbed  lattice.  It 
remains  to  be  demonstrated  that  such  Wannier  functions  exist.  The  Wannier 
functions  as  we  introduced  them  (see  equation  (4)) are  defined  as  Fourier  coeffi 
clents  of  the  Bloch  functions  which  are  periodic  In  k-space.  The  quasi¬ 
momentum  Is  defined  by  the  translation  operator  [15]  only  for  periodic 
potentials.  At  a  heterojunction,  the  periodicity  Is  broken  and  no 
quasi -momentum,  Bloch  functions  or  Wannier  functions  can  be  defined.  But  one 
can  attempt  to  generalize  the  entitles,  as  the  periodicity  is  essentially 
broken  only  at  the  junction.  The  most  obvious  approach  is  to  generalize  the 
Wannler  functions  which  are  strongly  localized  on  each  lattice  site  and  there¬ 
fore  will  be  strongly  perturbed  only  at  the  junction  itself. 

The  existence  of  generalized  Wannler  functions  (6WF)  has  been  theoreti¬ 
cally  established  by  Kohn  and  Onffroy  [21]  for  one-dimensional  structures. 
These  Wannler  functions  can  be  labeled  by  their  generalized  band  index;  they 
account  however  for  the  presence  of  interface  states  in  the  forbidden  band  gap 
[22].  We  shall  denote  these  generalized  Wannier  functions  a(n,x).  They 
span  the  same  function  space  as  the  eigenfunctions  of  the  perturbed  lattice 
for  the  same  generalized  band. 
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The  Index  n  will  be  taken  here  as  the  site  Index,  although  depending  on 
the  potential  well  there  may  be  more  than  one  or  even  no  GWF's  associated  with 
the  Interface  sites.  The  GWF's  are  localized  about  the  lattice  site  (or  the 
periodic  potential  minima).  Rehr  and  Kohn  have  shown  [23]  that  they  can  be 
exponentially  localized.  An  essential  property  of  the  GWF's  Is  that  they  form 
an  orthonormal  set 

<a(n) |a(m)>  «  /  a(n,x)  a(m,x)  dx  »  «n  m 

The  eigenfunctions  are  linear  combinations  of  a(n,x) 

T(x,E)  -  l  f(n,E)  a(n,x) 
n 

with  f(n,E)  the  solution  of 

E  f(n,E)  -  l  H(n,m)  f(m,E)  (12) 

m 

where  H(n,m)  Is  the  Hamiltonian  matrix  element  given  by 
H(n,m)  »  <a(n) |H|a(m)>. 

Here  H  Is  the  total  Hamiltonian  of  the  heterojunction  given  by 
H  -  H0  +  U(x) 

with  Hq  the  unperturbed  crystal  potential  and  U(x)  the  heterojunction 
effective  potential. 

Equation  (12)  Is  similar  to  equation  (7)  except  that  we  use  the  GWF 
Instead  of  the  unperturbed  Wannier  functions.  The  knowledge  of  the  GWF  weight 
f(n,E)  is  sufficient  in  most  applications  and  equation  (12)  can  be  solved  if 
the  matrix  element  H(n,m)  is  known. 


The  derivation  of  the  matrix  element  Is  the  critical  part  of  this  method. 
One  approach  consists  of  first  deriving  the  GUF's  at  each  site.  An  Important 
property  of  the  GWF's  demonstrated  by  Kohn  and  Onffroy  states  that  for  lattice 
sites  far  from  the  Interface  the  GWF's  approach  the  unperturbed  WF 
exponentially  In  n, 

a(n,x)  ■  w(x-na)  for  |n|  >>  0 

assuming  the  interface  at  n*0.  Only  a  few  GWF's  are  found  to  be  perturbed  or 
to  differ  from  the  WF  as  was  reported  by  Gay  and  Smith  [22],  It  follows  that 
a  few  lattice  parameters  away  from  the  Interface  the  Wannier  picture  developed. 
In  Sect1on3.1  holds  and  the  only  remaining  unknown  terms  are  the  interface 
Hamiltonian  matrix  terms. 

A  variational  method  for  the  derivation  of  the  GWF  has  been  proposed  by 
Kohn  and  Onffroy  for  one-dimensional  potentials  and  was  tested  by  [22].  A 
three-dimensional  treatment  Is  required  for  the  heterostructures  of  Interest 
here.  The  derivation  of  the  Wannier  functions  In  three  dimensions  being  quite 
Involved  [13],  the  derivation  of  the  more  complicated  GWF  will  not  be 
considered  as  a  practical  approach. 

For  our  present  purposes,  we  shall  use  for  the  Hamiltonian  matrix  simple 
approximations  derived  from  the  general  properties  of  the  Wannier  Functions 
discussed  in  Sections  3.1  and  3.2.. 

The  essential  result  brought  out  by  the  GWF  formalism  is  that  as  long  as 
a  generalized  band  of  energy  can  be  assumed.  It  is  possible  to  lump  all  the 
bandmixing  effects  Into  the  Interface  Hamiltonian  matrix  elements,  whereas  the 
wavefunctlon  f(n)  remains  given  by  an  equivalent  Hamiltonian  recurrence 
equation. 


3.3  SIMPLE  EXAMPLES 

We  shall  now  consider  two  simple  applications  of  the  Wannler picture. 

The  first  consists  of  the  evaluation  of  the  resonant  states  In  a  square 
quantum  well  .5  eV  deep  using  GaAs  bandstructure.  One  expects  to  observe  the 
Interaction  of  the  upper  valley  and  central  valley  states  of  GaAs  for  such  a 
deep  well  [24]. 

These  bandstructure  data,  kindly  furnished  by  Professor  Karl  Hess  of  the 
University  of  Illinois,  were  given  at  a  set  of  sampled  points  In  the  Brillouin 
zone.  Using  the  sampling  theorem  we  reconstituted  the  bandstructure.  The 
bandstructure  was  In  addition  fitted  at  r  to  .068 mQ.  This  consequently 
slightly  altered  the  smoothness  of  the  upper  valley.  The  resulting  band- 
structure  of  GaAs  along  a  Is  shown  In  Figure  1.  The  amplitudes  of  the  Fourier 
coefficients  of  the  GaAs  bandstructure  along  a  are  plotted  in  Figure  2.  We 
present  In  Figure  3  the  locus  of  the  amplitude  of  the  roots  exp  [jki(E)a)] 
that  solve  the  equation  E(ki)  *  E  where  E(k)  is  the  bandstructure  along  a. 

One  obtains  twenty  roots.  The  roots  1  and  2  correspond  to  Bloch  wave 
functions  of  the  central  valley  at  r.  The  roots  3  and  4  correspond  to  the 
Bloch  waves  of  the  upper  valley  centered  on  X.  The  unexpected  behavior  at  X 
of  the  roots  5,  6,  7  and  8  is  due  to  the  fact  that  the  upper  valley  minimum 
does  not  occur  exactly  at  X.  The  shift  of  the  X  minimum  is  an  artifact  of  the 
spectrum  truncation.  Since  the  bandstructure  is  slightly  distorted,  we  do  not 
expect  these  roots  to  contribute  much.  There  are  finally  twelve  remaining 
roots  with  a  very  high  damping  rate.  These  roots  contribute  only  within  the 
matching  area.  One  observes  in  Figure  3  the  evanescent  Bloch  waves  being 
transformed  into  propagating  Bloch  waves  at  r(0  eV)  and  X  (.3  eV). 


The  Wannler  picture  calls  for  the  matrix  element  of  the  quantum  well 
potential.  For  trial  purposes  we  evaluated  numerically  (with  the  Gauss 
Legendre  Integration  technique)  Its  matrix  element  using  the  WF's  the  weakly 
localized  WF's  of  the  empty  lattice 

-(X-/P)  ■  (f)1/a 

It  Is  most  enlightening  to  consider  the  diagonal  terms  plotted  In  Figure  4 
(with  an  arbitrary  energy  scale).  One  observes  that  these  diagonal  terms 
approximately  follow  the  assumed  square  potential  even  for  such  a  narrow  well. 
The  energy  difference  (a  few  percent)  Is  distributed  In  the  non-diagonal  terms 
and  contributes  somewhat  to  the  matching.  Using  the  complete  matrix  element 
we  found  that  sufficiently  accurate  results  could  be  generated  by  directly 
reducing  the  matrix  element  to  diagonal  terms  corresponding  to  the  sampling  of 
the  assumed  potential  at  the  lattice  sites  (see  Equation  9). 

We  have  developed  an  algorithm  to  solve  the  eigenvalue  problem  and 
evaluated  the  locus  of  the  resonant  energy  level  of  the  square  quantum  well  of 
depth  .5  eV,  as  a  function  of  the  well  width  L.  The  result  is  plotted  in 
Figure  5  (circles)  and  compared  with  the  energy  level  locus  (full  line)  of  the 
effective  mass  approximation. 

It  should  be  pointed  out  that  in  our  picture,  the  well  width  can  only  be 
varied  by  increments  of  a  half-lattice  parameter  for  theflOOldi rection. 

We  label  the  first  resonant  state  1.  For  a  wide  potential  well,  both  the 
WF  and  the  continuum  methods  yield  the  same  results,  since  the  resonant  level 
lies  at  low  energy  where  the  effective  mass  is  quite  accurate.  For  a  narrower 
well,  one  observes  a  divergence  of  the  two  energy  plots  as  the  effective  mass 
approximation  becomes  Inadequate  at  higher  energy. 


Ill 


GaAs  QUANTUM  WELL  (Wed  depth;  .5  eV) 

FIGURE  5  Energy  levels  of  a  GaAs  quantum  well,  as  a  function 
of  the  well  width. 
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For  a  well  smaller  than  five  lattice  parameters,  the  first  resonant  state 
occurs  above  .3  eV.  This  resonant  state  results,  then,  from  a  superposition 
of  the  central  valley  Bloch  wave  and  the  upper  valley  Bloch  wave.  Indeed,  one 
observes  the  interference  of  the  energy  locus  (labeled  a)  of  the  first 
resonant  level  of  the  upper  valley  with  the  energy  locus  1  of  the  central 
valley.  This  destroys  the  smooth  variation  of  these  energy  loci  as  they 
interact.  The  same  phenomenon  occurs  also  as  the  second  and  third  resonant 
l^vel  loci  (denoted  respectively  2  and  3)  cross  the  X  relative  minimum  (.3 
sV),  A  more  thorough  discussion  of  this  interference  effect  has  been  recently 
reported  by  Yla-Chung  Chang  and  D.Z.-Y.  Ting  [24], 

The  phenomenon  just  described  constitutes  an  example  of  what  may  be 
called  a  full  bandstructure  effect  which  can  be  handled  by  the  generalized 
Wannier  picture.  The  algorithm  developed  can  be  used  to  study  any  type  of 
variation  of  the  bandgap.  As  a  second  example,  we  consider  the  simple  case  of 
a  superlattice  made  of  square  wells  such  as  In  the  Kronig-Penney  proglem  [25]. 
We  assume  a  simple  tight-binding  bandstructure  E(k)  =  A[1  -  cos  lea]  with  the 
effective  mass  at  k*0  selected  to  be  .12  mo.  The  well  depth  is  .2  eV,  the 
width  and  separation  of  the  wells  is  chosen  to  be  10a  so  that  the  superlattice 
parameter  is  20a.  This  leads  to  two  resonant  levels  in  an  individual  well. 

In  Figure  6  we  present  the  result  for,  successively,  one,  two,  and  twenty 
wells.  The  discrete  bandstructure  of  the  20-well  superlattice  is  plotted  in 
Figure  7  (circles).  This  is  to  be  compared  to  the  bandstructure  (full  line) 
of  the  Kronig-Penney  model  in  the  effective  mass  approximation. 

The  low-lying  bands  (band  1)  have  a  similar  shape  but  are  separated  by  a 
shift  of  1  meV.  This  small  shift  originates  from  the  use  of  the  approximate 
diagonal  matrix  element  for  the  superlattice  potential. 


FIGURE  6  Superlattice  formation 
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FIGURE  7  Superlattice  bandstructure  (o:tight  binding  super 
lattice  of  20  quantum  wells);  (full  line:  Kronig 
Penney  model.) 


The  upper  bands  (band  2),  however,  differ  appreciably.  Part  of  the 
difference  Is  due  to  the  deviation  from  parabolicity  of  the  cosine  band 
and  results  in  a  much  higher  conduction  band  maximum  for  the  effective  mass 
approximation.  The  remaining  difference  Is  essentially  due  to  the  finite 
character  of  the  superlattice  composed  of  only  twenty  wells. 

More  sophisticated  structures  can  be  studied  with  the  algorithm 
developed,  but  this  was  not  our  purpose,  as  we  Intended  to  test  the 
generalized  Wannler  picture  using  simple  examples  for  which  closed-form 
solutions  are  available. 

3.4  "TRUE"  HETEROSTRUCTURES 

The  Wannler  picture  developed  above  and  Its  generalization  does  not  apply 
so  far  to  bandstructures  varying  with  space.  Let  us  point  out  right  away  that 
the  Idea  of  spatially  varying  bandstructures  Is  an  111-deflned  concept,  which 
can  be  made  more  rigorous  within  the  framework  of  the  generalized  Wannler 

model. 

A  one-dimensional  approach  will  breakdown  for  the  proper  description  of  a 
spatially  varying  bandstructure,  since  the  variation  In  the  bandstructure 
occurs  In  the  three-dimensional  k-space  and  therefore  requires  a 
three-dimensional  matching  of  the  Wannler  envelope. 

We  shall,  however,  for  simplicity  partially  retain  the  one-dimensional 
picture  by  assuming  the  spatially  varying  bandstructure  to  be  given  by 

ki2 

E(ic,x)  «  E(k,x)  ♦  Yn i*“ 

where  m*  Is  a  transverse  effective  mass  which  is  not  varying  In  space.  Such 


a  bandstructure  does  not  correspond  to  any  realizable  Interface  but 
constitutes  a  convenient  model.  In  this  model  vie  use  the  one-dimensional  GWF 
formalism  to  solve  for  the  eigenstates  of  this  structure.  A  set  of  orthogonal 
GWF's  noted  a(nfx)  can  be  defined  as  In  Section 3.2.  These  GUF's  can  be 
exponentially  localized  at  each  lattice  site  n.  For  the  case  of  an  abrupt 
heterojunction,  the  GUF's  a(n,x)  approach  exponentially  with  n  the  UF's  of 
semiconductors  1  and  2,  for  lattice  sites  far  from  the  Interface  (at  n*0). 

A  few  lattice  sites  away  from  the  Interface,  inside  semiconductor  1,  a  GWF 
closely  resembles  the  WF  of  the  unperturbed  semiconductor  1;  the  orthogonality 
with  the  GUF's  of  semiconductor  2  across  the  interface  essentially  affects 
only  the  exponentially  decaying  tall  of  this  GWF.  Finally,  one  Is  again  led 
to  the  generalized  one  band  formalism  of  Section 3.2 with  the  GWF  envelope 
given  by  the  one-dimensional  Wannler  recurrence  equation  (7). 

The  derivation  of  the  Hamiltonian  matrix  elements  H(n,m)  Is  again 
critical  to  the  analysis.  We  favor  a  direct  evaluation  of  these  elements 
which  remains  to  be  devised.  In  Its  absence  alternative  approaches  can  be 
conceived,  and  we  shall  demonstrate  such  an  approach  for  a  simple  system. 

First  consider  the  trial  Hamiltonian 

H(n,m)  »  E(n-m,n)  +  U(n,m) 

where  E(n-m,n)  is  the  bandstructure  Fourier  coefficient  which  is  assumed  to 
vary  In  space  with  the  lattice  site  n.  Such  a  Hamiltonian  which  might  be 
proposed  for  a  smoothly  spatially-varying  bandstructure  is  not  correct  as 
the  Hamiltonian  is  no  longer  hermitean: 


This  arises  as  a  consequence  of  using  the  Intrinsically  Ill-defined  concept  of 
spatially  dependent  bandstructures.  The  GWF  formalism  Is  capable  of 
accounting  for  the  breakdown  of  the  bandstfucture.  The  hermltlclty  Is  Indeed 
compatible  with 

E(n-n,n)  #  E(n-m,n) 

It  follows  that  the  related  bandstructure  at  the  lattice  site  n 


E(k.n)  -  l  E(l,n)  e 


Is  now  Imaginary  and  therefore,  as  It  should,  has  lost  Its  usual  physical 
meanl ng. 

For  an  abrupt  heterojunction,  the. use  of  GWF  formalism  Is  reduced  to  the 
derivation  of  matching  rules.  Indeed  away  from  the  Interface  the  Wannler 
picture  of  Section  3.1  holds  and  the  only  unknown  terms  are  the  Interface 
Hamiltonian  element  H(n,m). 

Consider  the  tight-binding  bandstructures  given  by 
Ei(k,x)  •  A  -  A  cos  ka 
E2(k,x)  ■  B  -  B  cos  ka 

Where  the  amplitude  A  and  B  are  respectively  related  to  their  effective  mass 


a  1,2 
*  mia2 


""  8  *  ^ 


For  simplicity  we  evaluate  the  matrix  element  H(n,m)  using  the  unperturbed 
WF's  on  both  sides  of  the  junction.  In  the  tight  binding  approximation  the 
only  remaining  unknown  Is  then  the  matrix  element  H(1,1+l)  »  -  C/2  assuming 
the  Interface  Is  located  somewhere  between  the  lattice  site  1  and  1+1.  The 
resulting  Hamiltonian  is  shown  around  the  Interface  in  Figure  8  for  U(n,m)»0. 


The  evaluation  of  C  depends  now  on  the  particular  matching  theory  on 

which  we  choose  to  rely.  There  has  been  much  published  recently  on  the 

matching  of  wavefunctlons  across  a  heterojunction.  The  GWF  picture  proposes  a 

general  method  which  enables  us  to  evaluate  the  wavefunctlon  In 

heterostructures,  but  calls  for  the  knowledge  of  H(n,m).  In  its  absence  we 

Intend  first  to  relate  this  method  to  other  matching  theories. 

Popular  techniques  are  the  effective-mass  matching  (see  [26]  and  the 

discussion  by  Kroemer  and  Ql-Gao  Zhu  [27]-[28]),  matching  developed  from  the 
♦  ♦ 

(k-p)  Kane  model  by  Bastard  [29],  White,  Margues  and  Sham  [30]-[31],  and  from 
tight-binding  arguments  by  White,  et  al.  [30]  and  Ando  and  Mori  [32].  A 
review  of  the  different  methods  Is  briefly  given  In  White  et  al.  [30].  We 
Intend  to  give  a  more  complete  discussion  of  the  matching  problem  In  a  forth¬ 
coming  paper. 

It  follows  from  Zhu  and  Kroemer's  work  [28]  that  the  effective-mass 
matching  Is  related  for  small  energies  to  the  geometrical  average  for  the 
"Ideal"  heterostructure  (type  I  [28])  considered  here 
C  -  (A-B)1/2 

From  Bastard's  work  [29],  we  have  C  given  by  the  average 
C  *  (A  +  B)/2. 

Both  expressions  actually  lead  to  the  same  results  when  the  variation  of 
effective-mass  Is  small  across  the  heterojunction.  In  Figure  9  we  compare  the 
effective  mass  matching  with  the  geometrical  average  matching  for  a  quantum 
well.  This  quantum  well  Is  made  of  a  layer  of  Ga.47In.53As  sandwiched  between 
Al.48ln.52As.  The  effective  masses  are  respectively  taken  as  .041  and  .075  m0 
for  GalnAs  and  AlInAs.  These  effective  masses  should  be  corrected  for  narrow 
wells  as  the  non-parabolicity  is  no  longer  negligible  (see  Welch,  Wicks  and 


-Ill- 


3:(  c\ 


100 


QM/AfeiAa  QUANTUM  WELL  (DEPTH;  3  «V) 


FIGURE  9  Energy  locus  as  a  function  of  well  width  for 

Ga  47ln  53AsZAl.48ln.5As  quantum  well  (0: results 
using  Wannier  algorithm  with  the  geometrical  average 
matching);  (full  line:  results  using  the  effective 
mass  matching). 


-112- 


Eastman  [33]).  The  discrete  energy  locus  is  seen  to  agree  with  the  continuous 
energy  locus  based  on  the  effective  mass  formalism. 

This  result  was  established  within  the  tight-binding  approximation.  A 
more  refined  scheme  will  extend  the  connection  rules  to  Include  higher  order 
overlaps. 

3.5  CONCLUSIONS 

We  have  introduced  a  generalized  one-dimensional  Wannier  picture.  This 
formalism  enables  us  to  handle  sharp  variations  in  potential  as  long  as  the 
potential  strength  Is  sufficiently  moderate  for  the  concept  of  a  generalized 
band  to  hold.  The  associated  band-mixing  effects  are  then  lumped  into  a 
generalized  one-band  Hamiltonian  matrix.  In  the  absence  of  knowledge  of  the 
GWF  one  can  often  rely  in  practice  on  the  direct  sampling  of  the  potential  at 
each  lattice  site. 

The  generalized  one  band  picture  applies  directly  to  modulation-doped 
structures  and  can  be  used  as  a  model  of  heterostructures  by  treating  the 
bandgap  variation  as  an  effective  potential.  Another  application  not  covered 
here  and  for  which  this  model  was  in  fact  developed  concerned  the  study  of 
devices  such  as  the  Zener  superlattice  oscillator  for  which  the  periodicity  in 
k-space  of  the  subband  together  with  Its  spatial  variation  is  fundamental  to 
the  device  operation  [5]. 

Finally,  practical  derivation  of  the  energy  levels  and  eigenfunctions  can 
be  performed  for  arbitrary  potentials  or  heterostructure  configurations  using 
the  same  algorithm. 

The  one-dimensional  generalized  Wannier  formalism  can  be  used  to  model 
true  heterostructures,  although  a  rigorous  treatment  requires  a  three- 
dimensional  GWF  analysis.  In  the  one-dimensional  picture  the  matching 
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of  the  wavefunction  across  the  Interface  Is  determined  by  the  Interface 
Hamiltonian  matrix  elements. 

The  derivation  of  Wannler  Functions  Is  In  general  a  complex  problem  due 
to  their  orthogonality  and  the  fact  that  they  are  not  eigenfunctions  of  the 
Hamiltonian.  The  best  procedure  would  be  to  evaluate  the  generalized  Fourier 
coefficients  of  the  generalized  bandstructure  directly  just  as  the 
bandstructure  Is  evaluated  directly  In  pseudopotential  calculations  using 
variational  techniques. 

An  alternative  approach  Is  to  derive  the  Interface  Hamiltonian  matrix 
directly  from  the  enforcement  of  physical  constraints.  Such  an  approach  was 
undertaken  by  Kroemer  [28]  and  Ando  and  Mori  [32]  for  maximized  matching 
(Ideal  heterojunction)  using  current  continuity.  We  shall  report  In  a 
forthcoming  paper  on  the  extension  of  these  matching  theories  to  higher  energy 
range  (non-effective  mass  case)  and  for  the  inclusion  of  higher  order  overlap 
(non  tight-binding  cases)  for  both  type  I  and  II  semiconductors  (see  [28]). 

The  concept  of  generalized  bands  described  here  is  largely  Intuitive  and 
has  been  used  so  far  as  such.  The  Introduction  of  the  generalized  Wannier 
functions  brings  support  to  the  concept.  One  interesting  feature  of  the 
generalized  Wannier  picture  Is  to  account  consistently  for  the  concept  of 
spatially  varying  bandstructure.  The  generalized  Wannier  picture  appears  then 
as  an  intermediate  method  situated  between  the  direct  study  at  the  lattice 
level  (e.g. ,  LCOA  techniques)  and  the  techniques  based  on  the  effective  mass 
approximation. 

We  believe  that  useful  insights  can  be  drawn  from  this  ability  of  the  GWF 
formalism  to  handle  the  concept  of  bandstructure  through  the  use  of  the  band- 
structure  Fourier  coefficients  together  with  their  spatial  dependence  on  the 


lattice  sites.  The  bandstructure  appears  as  a  physical  entity  established  In  a 
few  lattice  parameters.  We  found  that  5  to  10  lattice  sites  were  sufficient 
for  a  close  representation  of  a  bandstructure  along  <100>.  This  Is  in  quali¬ 
tative  agreement  with  the  Wannler  functions  evaluated  by  Kane  and  Kane  [13]. 
Energy  bands  calculated  from  these  Wannler  functions  were  accurate  to  .1  eV 
when  the  Wannler  Function  overlaps  were  evaluated  to  twelfth  neighbors  (459 
individual  bonds).  The  spatial  establishment  of  a  small  effective  mass 
requires  a  larger  number  of  sites  since  a  large  number  of  Wannler  function 
overlaps  or  bandstructure  Fourier  coefficients  Is  required  to  fit  a  small 
effective  mass  while  accurately  reproducing  the  remainder  of  the  band- 
structure. 
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4.  LONG  RANGE  ORDERING 


For  a  semiconductor  alloy,  A  B,  C,  regular  solution 

A  X 

theory  indicates  that  the  excess  enthalpy  or  heat  of  mixing 
is  given  by 

a  Hm  =  x(l-x) n, 

where  Q  is  an  interaction  parameter  which  contains  the  physics 
of  the  mixing  process.  For  positive  values  of  Q  we  expect 
to  have  cation  clustering  above  and  phase  separation  below 
a  critical  temperature, 

^c  k41n2  * 

For  negative  values  of  a  we  expect  a  tendency  toward  cation 
ordering  above  and  compound  formation  below  T  . 

4.1  RECENT  ALLOY  FORMATION  MODELS 

There  have  recently  been  several  calculations  of  Q  for 
III-V  alloys  based  on  the  idea  that  most  or  all  of  the  heat 
of  mixing  is  due  to  bond  distortion.  Fedders  and  Muller  [1] 
calculate  n  under  the  assumption  that  it  is  due  to  the  changes 
in  the  bond  lengths .  a  and  b,  of  the  binary  compounds,  AC  and 
BC,  required  to  form  a  virtual  crystal  lattice.  This  gives 


n  =  x 


i 


where  X  can  be  determined  from  the  elastic  properties  of  the 
material.  Mikkelsen  and  Boyce  [2]  have  found  from  EXAFS 


measurements,  however,  that  the  AC  and  BC  bonds  in  alloys 
are  within  20  percent  of  their  values,  a  and  b,  in  the  binary 
compounds ,  which  reduces  X  by  about  20. 

To  maintain  essentially  binary  bond  lengths,  a  and  b, 
in  a  virtual  crystal  alloy  the  bond  ancles  must  also  be  distorted 
so  that  the  factor  of  20  reduction  above  is  an  under  estimate 
of  A.  Mikkelson  [3]  has  taken  into  account  both  bond  length 
and  angle  distortions  to  obtain  values  of  ft  in  good  agreement 
with  experimental  values  of  ft  from  pseudo-binary  phase  diagrams. 
Although  Muller  [4]  has  shown  that  bond  distortion  should 
reduce  disorder  below  that  expected  for  a  random  cation  distribu¬ 
tion  of  A  and  B,  all  elastic  models  necessarily  result  in 
a  positive  heat  of  mixing  since  strain  energy  is  positive 
definite. 

Zunger  et  al.  [5]  have  included  charge  transfer  effects 
with  bond  distortion  in  their  model  of  semiconductor  alloy 
formation.  They  find  that,  although  the  negative  excess 
energy  from  charge  transfer  is  small  in  comparison  to  the 
positive  excess  heat  of  mixing  from  bond  distortion,  the 
chemical  charge  transfer  provides  a  means  of  stabilizing 
long  range  order  and  compound  formation.  This  is  because 
ordered  alloys  or  compounds  can  accomodate  bond  distortion 
with  less  strain  energy  than  disordered  alloys  with  large 
strain  energy  associated  with  random  cation  distribution. 
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As  an  example,  for  the  alloy  InQ  5GaQ  5P  they  obtain 
3  =  +0 . lleV/ atom-pair  which  is  the  same  as  obtained  by  Mikkelsen. 
For  the  compound  InGaP2  Zunger  et  al.  obtain  an  equivalent 
excess  energy  parameter  =  -  0 . 047eV/atom-pair .  On  this  basis 
we  can  qualitatively  compare  alloys  with  In  and  Ga  cations  to 
those  with  Al  and  Ga  cations.  That  is,  alloys  with  In  and  Ga 
cations  are  expected  to  have  relatively  large  bond  distortion 
(atomic  size  differences)  and  large  charge  transfer  (electro¬ 
negativity  differences)  compared  to  alloys  with  Al  and  Ga  cations. 
Thus,  the  small  positive  or  zero  heat  of  mixing  in  alloys  with 
Al  and  Ga  cations  is  proportionally  offset  by  small  negative 
or  zero  excess  energy  due  to  chemical  interaction. 

4 . 2  X-RAY  MEASUREMENTS 

As  mentioned  in  our  last  annual  report  [6]  we  have  obtained 
some  evidence  for  long  range  ordering  in  InxGai_xp  grown  on  GaAs 
by  liquid  phase  epitaxy.  This  is  indicated  in  Figure  1  where 
we  show  back-reflection  Laue  photographs  for  (a)  a  {211} Ga 
oriented  GaAs  substrate  and  (b)  an  InQ  4gGaQ  51P  lattice-matched 
epitaxial  layer  using  a  Cu  target.  In  these  measurements 
additional  spots  in  the  photograph  for  the  epitaxial  layer 
could  easily  be  interpreted  as  evidence  for  long  range  order. 

What  we  see  instead  are  enhanced  intensity  0  46  and  0  6  4  spots. 
These  enhanced  intensity  spots  were  also  observed  in  about 
eight  other  In  Ga,  P  layers  with  x  between  0.49  and  0.52,  but 
were  not  observed  in  one  A1Q  4Gao  gAs  laYer  we  measured.  When 
Laue  photographs  f or  an  InQ  4gGag  51P  layer  were  taken  with  a 
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Figure  1  Back-reflection  Laue  photographs  for  (a)  a  {211} Ga  orientated 
GaAs  substrate  and  (b)  an  Ino.49Gao.51P  lattice-matched 
epitaxial  layer.  The  primary  X-ray  beam  from  a  Cu  target 


Ho  target ,  the  0  4  6  and  0  6  4  spots  were  not  enhanced  but 
additional  210,.12,  2.10  12,,  2  12  10,  and  2  12,10,  were  observed 


Table  I.  Reflections  which  satisfy  the  Bragg  condition  for 
characteristic  Ka  lines  near  the  lattice  constant 
of  In_  ,Ga«  -P. 


Equivalent 

Cubic 

Reflection 


212  10 


4109 


Structure  Factors 


Sphalerite  Chalcopyrite  Layered 


f In”fGa 


fIn+fGa“2fp  fIn+fGa“2fp  fIn+fGa“2fp 


fIn”fGa 


fIn  +  fGa+2fp  fIn+fGa+2fp  fIn+fGa+2fp 


fIn”fGa 


fIn~fGa 


Characteristic 

Bragg 

Condition 


2nd  Order 
Cu  Ka 


1st  Order 
Cu  Ka 


2nd  Order 
Mo  Ka 


1st  Order 
Mo  Ka 


1st  Order 
Cu  Ka 


1st  Order 
Cu  Ka 


In  Table  I  we  compare  the  structure  factors  for  several 
reflections  in  the  sphalerite  structure  and  two  of  the  most 
likely  ordered  structures  for  InGaP2,  chalcopyrite  (space 
group  14  2d)  and  primitive  layered  (space  group  P4  2m)  structures. 
First,  note  that  superlattice  structure  factors  are  small 


Ga)  and  will  be  observed  only  when  they  satisfy  the 


Bragg  condition,  2dn^  sin  9  =  n  A ,  for  one  of  the  high  intensity 
characteristic  x-ray  lines.  For  example,  the  2  12  7  and  4 109 
(equivalent  cubic  notation)  chalcopyrite  reflections  satisfy 
the  Bragg  condition  in  1st  order  for  Cu  Ka  (A  =  1.5418A)  lines 
near  the  lattice  constant  of  GaAs  (5.6532A).  Since  we  easily 
observe  these  reflections  in  our  chalcopyrite  ZnSnP2  (5.6407A) 
epitaxially  grown  on  GaAs,  we  conclude  that  the  In  Ga,  P  is 

a  X  A 

not  ordering  in  a  chalcopyrite  structure. 

Second,  notice  that  the  064  set  of  reflections  are  allowed 
in  all  three  crystal  structures  and  satisfy  the  Bragg  condition 
in  1st  order  for  Cu  Ka.  If  the  In  Ga,  P  has  a  layered 

A  X  A 

structure,  these  064  reflections  could  be  enhanced  by  032 
reflections,  which  satisfy  the  Bragg  condition  in  2nd  order 
for  Cu  Ka.  In  a  similar  manner  the  2  12  10  reflections,  which 
satisfy  the  Bragg  condition  in  1st  order  for  Mo  Ka(A=  0.7107A) , 
could  be  enhanced  by  165  reflections,  which  satisfy  the  Bragg 
condition  in  2nd  order  for  Mo  Ka . 

Using  tabulated  atomic  scattering  factors  [ 7 ]  we  can 
calculate  the  intensity  of  the  064  Laue  spots  that  would 
be  expected  with  different  materials  and  crystal  structures. 
Taking  into  account  032,  064,  096,  and  0128  reflections, 
the  following  reflected  intensity  ratios  are  obtained. 

InGaP2  InQ  ^Ga0  5?  InGaP2 

InQ  5GaQ  5P  GaAs  GaAs 

3.5  60  220 

Although  this  calculation  neglects  absorption  and  secondary 

extinction,  it  shows  that  the  064  enhancement  shown  in 


125- 


Figure  1  could  be  simply  due  to  the  difference  between 
sphalerite  InQ  5GaQ  5P  and  GaAs  and  not  layered  InGaP2 
and  GaAs. 

4 . 3  LAYERED  STRUCTURE 

Recently,  Kuan,  Kuech  and  co-workers [ 8 ]  have  reported 
long-range  order  in  AlxGa1-xAs.  They  claim  to  see  the  most 
pronounced  evidence  for  ordering  from  A1Q  75GaQ  25As  9rown 
on  {110} GaAs  at  700°C  by  organo-metallic  vapor  phase  epitaxy. 
Electron  diffraction  measurements,  which  are  sensitive  to 
surface  structure,  were  used  to  indicate  the  presence  of 
a  layered  crystal  structure  ( the  P  4  2  m  space  group  discussed 
above)  with  c-axes  orientations  in  the  growth  plane.  No 
significantly  different  physical  properties,  however,  were 
obtained  for  this  material.  Although  this  report  is  probably 
premature  for  the  reasons  indicated,  there  is  some  merit 
to  their  claim. 

A  fairly  definitive  x-ray  experiment  which  would  indicate 
the  presence  of  a  layered  crystal  structure  is  shown  schematic¬ 
ally  in  Figure  2.  The  geometry  of  the  experiment  is  the 
same  as  that  used  to  obtain  the  back-reflection  Laue  photographs 
shown  in  Figure  1.  The  layered  face  centered  cubic  cation 
sublattice  with  cations  as  open  and  closed  circles  and  [001] 
c-axes  is  shown  as  observed  in  the  [100]  direction.  Figure  2  (a) 
shows  that  for  this  c-axis  orientation  reflections  from  the 
(046)  planes  between  the  (023)  planes  do  not  cancel  out  (023) 
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reflections.  Thus,  both  0  4  6  and  02  3  superimpose  on  the 
bottom  of  the  film.  Figure  2  (b)  shows,  however,  that  for 
the  same  c-axis  orientation  reflections  from  the  (064)  planes 
cancel  out  reflections  from  the  (032)  planes.  Thus,  only 
0  6  4  is  observed  on  the  top  of  the  film.  A  significant 
difference  in  intensity  of  these  spots  on  the  film  (we  calcu¬ 
late  3.5  for  In  and  Ga  cations),  therefore,  would  indicate 
a  layered  structure  with  one  c-axis. 

If  the  c-axis  in  Figure  2  were  in  the  (010]  direction, 
one  can  see  that  the  '0  6  4  and  0  3  2  at  the  top  of  the  film 
would  superimpose  while  the  023  reflection  would  be  cancelled 
by  064  at  the  bottom.  For  the  c-axis  in  the  [100]  direction 
both  0  3  2  and  0  2  3  reflections  would  be  cancelled,  since 
the  layered  structure  looks  like  a  sphalerite  structure  to 
an  [Oil]  x-ray  beam  orthogonal  to  the  c-axis.  The  examination 
of  a  tetragonal  structure  with  fixed  primary  x-ray  beam 
orientation  and  three  c-axis  orientations  is  equivalent  to 
the  examination  of  a  tetragonal  structure  with  fixed  c-axis 
orientation  and  three  primary  x-ray  beam  orientations.  The 
latter,  of  course,  is  experimentally  achievable  and  we 
summarize  potential  observations  in  Table  II . 


Beam 

Directions 

Enhanced 

0  6  4  Spots 

P  4  2m 
c-axes 

[Oil] 

0  4  6  — 

[101] 

_  4  0  6 

[001] 

[110] 

—  — 

[Oil] 

0  4  6 

0  6  4 

[001] 

[101] 

—  4  06 

and 

[110] 

4  6  0  — 

[010] 

[Oil] 

0  46 

064 

[001] 

[101] 

.... 

[010] 

6  0  4 

4  0  6 

and 

[110] 

6  4  0 

4  6  0 

[100] 

The  salient  feature  of  the  experiments  indicated  in 
Table  II  is  that,  to  determine  the  number  of  c-axes  in  a 
layered  sample,  x-ray  measurements  must  be  taken  with  the 
primary  beam  in  at  least  two  of  the  <  110  >  directions 
indicated.  At  the  present  time  we  have  looked  at  our 
InJCGa1_xP  samples  with  the  primary  beam  in  only  one  <110> 
direction.  For  these  samples  no  asymmetry  in  064  spot 
intensity  has  been  observed  which  could  not  be  attributed 


to  slight  misorientation  of  the  primary  beam. 


4 . 4  CONCLUSIONS 

From  the  discussion  above  we  conclude  that  charge  transfer 
effects  in  III-V  alloys  may  provide  a  mechanism  for  stabilizing 
equilibrium  long  range  order  and  compound  formation.  The 
resulting  total  negative  excess  energy  involved,  however, 
appears  to  be  sufficiently  small  that  low  growth  temperatures 
and  kinetic  effects,  such  as  substrate  orientation,  may  be 
critical  factors.  This  appears  to  be  the  case  for  the  recently 
reported  ordering  in  Al  Ga,  As  and  should  also  be  true 

X"a 

for  alloys  with  larger  atomic  size  and  electronegativity 
differences. 
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